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The direct photon spectra and flow (v 2 , vf) in heavy-ion collisions at SPS, RHIC and LHC energies 
are investigated within a relativistic transport approach incorporating both hadronic and partonic 
phases - the Parton-Hadron-String Dynamics (PHSD). In the present work, four extensions are 
introduced compared to our previous calculations: (i) going beyond the soft-photon approximation 
(SPA) in the calculation of the bremsstrahlung processes meson + meson —>■ meson + meson + 7 , 

(ii) quantifying the suppression due to the Landau-Pomeranchuk-Migdal (LPM) coherence effect, 

(iii) adding the additional channels V + N ^ N + j and A ^ N + y and (iv) providing PHSD 
calculations for Pb-fPb collisions at ^^/SNN = 2.76 TeV. The first issue extends the applicability of the 
bremsstrahlung calculations to higher photon energies in order to understand the relevant sources in 
the region pr = 0.5—1.5 GeV, while the LPM correction turns out to be important for pr < 0.4 GeV 
in the partonic phase. The results suggest that a large elliptic flow V 2 of the direct photons signals 
a significant contribution of photons produced in interactions of secondary mesons and baryons in 
the late (hadronic) stage of the heavy-ion collision. In order to further differentiate the origin of the 
direct photon azimuthal asymmetry (late hadron interactions vs electromagnetic fields in the initial 
stage), we provide predictions for the photon spectra, elliptic flow and triangular flow V 3 {pt) of 
direct photons at different centralities to be tested by the experimental measurements at the LHG 
energies. Additionally, we illustrate the magnitude of the photon production in the partonic and 
hadronic phases as functions of time and local energy density. Finally, the ’cocktail’ method for an 
estimation of the background photon elliptic flow, which is widely used in the experimental works, 
is supported by the calculations within the PHSD transport approach. 

PACS numbers: 25.75.-q, 13.60.Le, 14.40.Lb, 14.65.Dw 


I. INTRODUCTION 

Direct photons are a powerful probe of the quark- 
gluon plasma (QGP) as created in ultra-relativistic nu¬ 
clear collisions. The photons interact only electro- 
magnetically and thus escape to the detector almost 
undistorted through the dense and strongly-interacting 
medium. Thus the photon transverse-momentum spec¬ 
tra and their azimuthal asymmetry carry information on 
the properties of the matter under extreme conditions, 
existing in the first few fm/c of the collisional evolution. 

On the other hand, the measured photons provide a 
time-integrated picture of the heavy-ion collision dynam¬ 
ics and are emitted from every moving charge - partons 
or hadrons. Therefore, a multitude of photon sources has 
to be differentiated in order to access the signal of inter¬ 
est. The dominant contributions to the inclusive photon 
production are the decays of mesons, dominantly pions, 
rj- and w-mesons. Experimental collaborations subtract 
the “decay photons” from the inclusive photon spectrum 
using a cocktail calculation [I|, Q and obtain the “direct” 
photons. 

In particular the direct photons at transverse momenta 
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(pt < 3 GeV) are expected to be dominated by ’’ther¬ 
mal” sources, i.e. the radiation from the strongly in¬ 
teracting Quark-Gluon-Plasma (sQGP) Q and the sec¬ 
ondary meson-bmeson and meson-l-baryon interactions in 
the hadronic phase 0,0. These partonic and hadronic 
channels have been studied within PHSD in detail in 
Refs. 0,0 at Relativistic-Heavy-Ion-Collider (RHIC) en¬ 
ergies. It was found that the partonic channels constitute 
up to half of the observed direct photon spectrum for very 
central collisions. Other theoretical calculations also find 
a significant contribution of the photons produced in the 
QGP to the direct photon spectrum [§-[1^. 


The low-pT direct photons probe not only the temper¬ 
ature sain of the produced QGD-matter, but also 
its (transport) properties, for instance, the shear viscos¬ 
ity. Using the direct photon elliptic flow V 2 (a measure 
of the azimuthal asymmetry in the photon distribution) 
as a viscosimeter was first suggested by Dusling et al. in 
Ref. [l^; this idea was later supported by the calcula¬ 
tions in Refs. [l3, [HI, [H, Ell • It was also suggested that 
the photon spectra and t >2 are sensitive to the collective 
directed flow of the system 0113, to the equation of 
state 0E3, to the possible production of a Glasma 


to the rate of chemical equilibration in the QGP 


21 

24| and to the asymmetry induced by the strong magnetic 


field (flash) in the very early stage of the collision |2£ 


However, the recent observation by the PHENIX Col- 
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laboration [I| that the elliptic flow V 2 {pt) of ’direct pho¬ 
tons’ produced in minimum bias Au-|-Au collisions at 
y/SNN = 200 GeV is comparable to that of the produced 
pions was a surprise and in contrast to the theoretical 
expectations and predictions. Indeed, the photons pro¬ 
duced by partonic interactions in the quark-gluon plasma 
phase have not been expected to show considerable flow 
because they are dominated by the emission in the initial 
phase before the elliptic flow fully develops. 

In Refs. in we have applied the PHSD approach 
to photon production in Au-l-Au collisions at -^/snn = 
200 GeV and studied the transverse momentum spec¬ 
trum and the elliptic flow V 2 of photons from hadronic 
and partonic production channels. The microscopic de¬ 
scription of the full collision evolution is calculated in 
the covariant off-shell transport approach PHSD. The 
degrees-of-freedom in the partonic and hadronic phases 
are strongly interacting dynamical quasi-particles and 
off-shell hadrons, respectively. 

It was found in Refs. @,0 that the PHSD calculations 
reproduce the transverse momentum spectrum of direct 
photons as measured by the PHENIX Collaboration in 
Refs. [11,133. The centrality dependence of the thermal 
photon yield in PHSD was predicted to be ~ A^part '"’ith 
the exponent a = 1.5, which is in good agreement with 
the most recent measurement of a = 1.48 ± 0.08 ± 0.04 
by the PHENIX Collaboration [3l|. Also, the PHSD de¬ 
scribed the data on the elliptic flow of inclusive and direct 
photons at the top RHIC energy. The strong V 2 of direct 
photons - which is comparable to the hadronic V 2 - in 
PHSD was attributed to hadronic channels, i.e. to meson 
binary reactions which are not subtracted in the data. As 
sources for photon production, we have incorporated the 
interactions of off-shell quarks and gluons in the strongly 
interacting quark-gluon plasma (sQGP) {q + q ^ g + ^ 
and q{q)+g —t q{q)+"f), the decays of hadrons (tt — >■ 7-I-7, 
T] —>■ 7-I-7, w —^ 7r-|-7, 77' —^ p-1-7, (/> —t ?? + 7, oi — 1 77-1-7) 
as well as their interactions 7r-|-7r —>■ p + ^, p + ir —>■ 7r-|-7, 
meson-meson bremsstrahlung m m —> TO-l-m-l- 7 ), 
meson-baryon bremsstrahlung {m + In 

the present work we additionally incorporate the two-to- 
two vector meson-|-nucleon interactions (V ^ ^+p/n 
and V + n ^ ^+pln) and the decay of the A-resonance 
A —>■ Nj. 

The photon production via bremsstrahlung in meson- 
meson and meson-baryon elastic collisions was found to 
be a very important source to interpret the data on 
the direct photon spectra and elliptic flow simultane¬ 
ously i0- In the previous works i0 , we have been 
calculating the photon bremsstrahlung from all elastic 
meson-meson and meson-baryon scatterings mi + m 2 , 
which occur during the heavy-ion collisions (including 
mi = TT, 7 , K, K, K^, K*, K*, , g',ui, p,4>,ai). Eor the 
calculation of the bremsstrahlung cross sections we have 
been applying the soft-photon approximation (SPA). 
Therefore the resulting yield of the bremsstrahlung pho¬ 
tons depended on model assumptions such as i) the cross 
section for the meson-meson elastic scattering (we as¬ 


sumed 10 mb for all meson species), ii) incoherence of the 
individual scatterings and iii) the soft-photon approxima¬ 
tion (i.e. low photon energy and low ^/s of the collision). 
The adequacy of the SPA assumption has been studied 
before in Ref. 32| and a theoretical uncertainty of up to 
a factor of 2 was found. 

In view of the importance of questions that can be 
answered by direct photon measurements, we have im¬ 
proved our implementation of the photon production 
in the PHSD transport approach, in particular of the 
bremsstrahlung photon production in the meson-meson 
collisions. In the scope of the present work, we depart 
from the assumptions above by i) using a covariant one- 
boson-exchange (OBE) chiral model for the tt interac¬ 
tions, ii) investigating the suppression due to the coher¬ 
ence of the photon emission with long wavelength (LPM 
effect) and iii) deriving and implementing the “exact” 
OBE cross sections for photon bremsstrahlung thus de¬ 
parting from the SPA. 

The results of our calculations are compared to the 
data from the SPS and RHIG Gollaborations to check 
if the earlier conclusions and interpretations hold. Ad¬ 
ditionally, we provide calculations for the photon pro¬ 
duction in Pb -I- Pb collisions at the energy of ^/snn = 
2.76 TeV (spectra, elliptic V 2 and triangular V 3 flow of 
direct photons). The comparison of these calculations to 
the future data of the ALICE Collaboration will be of 
great interest since the preliminary data [1 ,113 indicate 
a significant direct photon signal at low px with a large 
elliptic flow. 


II. PHOTON PRODUCTION IN HEAVY-ION 
COLLISIONS WITHIN THE PHSD APPROACH 

The PHSD model [13, HH] is an off-shell transport ap¬ 
proach that consistently describes the full evolution of a 
relativistic heavy-ion collision from the initial hard scat¬ 
terings and string formation through the dynamical de¬ 
confinement phase transition to the quark-gluon plasma 
as well as hadronization and the subsequent interactions 
in the hadronic phase. The two-particle correlations re¬ 
sulting from the finite width of the parton spectral func¬ 
tions are taken into account dynamically in the PHSD by 
means of the generalized off-shell transport equations [13 
that go beyond the mean field or Boltzmann approx¬ 
imation [ 13 . The transport theoretical description of 
quarks and gluons in the PHSD is based on the Dynami¬ 
cal Quasi-Particle Model [13 (DQPM) for partons that is 
constructed to reproduce lattice QCD (IQCD) results for 
the entropy density, energy density and pressure as func¬ 
tions of temperature for the quark-gluon plasma in ther¬ 
modynamic equilibrium. In the hadronic sector, PHSD 
is equivalent to the Hadron-String Dynamics (HSD) ap¬ 
proach [ 39 I-I 41 II. For details about the DQPM model and 
the off-shell transport approach we refer the reader to the 
review Ref. [13 • 

We stress that a non-vanishing width in the partonic 
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spectral functions is the main difference between the 
DQPM and conventional quasiparticle models Its 
influence on the collision dynamics is essentially seen in 
the correlation functions. For instance, in the stationary 
limit, the correlation involving the off-diagonal elements 
of the energy-momentum tensor defines the shear vis¬ 
cosity 77 of the medium [i^. Here a sizeable width is 
mandatory to obtain a small ratio of the shear viscosity 
to entropy density 77 /s, which results in a roughly hydro- 
dynamical evolution of the partonic system in PHSD . 
The finite width leads to two-particle correlations, which 
are taken into account dynamically by means of the gen¬ 
eralized off-shell transport equations [s^, going beyond 
the mean field approximation [s^, ■ It has been shown 

in [ 4 ^ that the final width (the imaginary part of the self 
energy) is demanded by the causality constraint on the 
propagator in the theory of the strongly interacting par¬ 
ticles as soon as the interaction leads to a sizable dressing 
mass squared (real part of the self energy). 

In the past the PHSD approach has provided a con¬ 
sistent description of the bulk properties of heavy-ion 
collisions - rapidity spectra, transverse mass distribu¬ 
tions, azimuthal asymmetries of various particle species 

- from low Super-Proton-Synchrotron (SPS) up to the 

LHC energies [3111; 113 • hadronic sector, PHSD 

is equivalent to the Hadron-String Dynamics (HSD) ap¬ 
proach [ 3 ^ HO, Ell , in which the photon production at 
top SPS energies has been investigated before in Ref. [13 
with an emphasis on the role of meson-meson interac¬ 
tions. The PHSD approach was also successfully used for 
the analysis of penetrating probes, such as charm [4l.[5l| 
and dilepton production from hadronic an d p artonic 
sources at SPS, RHIC and LHC energies [^.1^. 

Indeed the calculations within the PHSD have re¬ 
produced the measured differential spectra of dileptons 
from heavy-ion collisions at SPS and RHIC energies (see 
Refs. El])- Also, it has been checked in Ref. [Sl] 
that the dilepton production from the QGP constituents 

- as incorporated in the PHSD [nm - agrees with the 
dilepton rate emitted by the thermalized QCD medium 
as calculated in the IQCD approach. Indeed, the decon- 
fined state of matter - created in heavy-ion collisions at 
RHIC (sl - fHll - was clearly seen in the dilepton yield 
above invariant masses of 1.2 GeV [sillsll. 

As sources of photon production - on top of the general 
dynamical evolution - we consider hadronic as well 
as partonic interactions. Let us first describe all 

the contributions, which consist of the photon produc¬ 
tion in the quark and gluon collisions, from the hadronic 
decays and the interactions of intermediate mesons pro¬ 
duced throughout the evolution of the nucleus-nucleus 
collision. 

1) Photons are radiated by quarks in the interaction 
with other quarks and gluons. In this sense, we differenti¬ 
ate two classes of processes: first the two-to-two reactions 


q + 9 + 1, 

q/q + 9 ^ q/q + 1- 


The implementation of the photon production by the 
quark and gluon interactions in the PHSD is based on 
the off-shell cross sections for the interaction of the mas¬ 
sive dynamical quasi-particles as described in li,ll 3 . In 
addition, photon production in the bremsstrahlung reac¬ 
tions q -I- q/g q + q/g -t- 7 is possible [^ . 

2) All colliding hadronic charges (meson, baryons) can 
also radiate photons by the bremsstrahlung processes: 


( 1 ) 

( 2 ) 

The processes ,0 have been studied within the PHSD 
in Refs. [ 1 , E3j while the processes ([2]) were added in 
Ref. [ 3 . The implementation of photon bremsstrahlung 
from hadronic reactions in transport approaches has 
been based until now on the ’soft photon’ approxima¬ 
tion (SPA). The soft-photon approximation |62l| relies on 
the assumption that the radiation from internal lines is 
negligible and the strong interaction vertex is on-shell 
and is valid only at low energy (and pr) of the produced 
photon. Since the relatively high transverse momenta of 
the direct photons {pr = 0.5 — 1.5 GeV) are most impor¬ 
tant for a potential understanding of the ” direct photon 
puzzle”, we extend here the accuracy of our calculations 
for the photons produced via the bremsstrahlung mech¬ 
anism beyond the applicability of the SPA. The details 
of the new model are given below in section IHI Al 

3) Additionally, the photons can be produced in binary 
meson meson and meson -I- baryon collisions. We con¬ 
sider here the direct photon production in the following 
2 —2 scattering processes: the meson-l-meson collisions 


TT -I- TT —p -I- 7, 

TT -I- p —>■ TT -I- 7, ( 3 ) 


and the meson-fbaryon collisions 


V + N 
where V = p, (j), w. 


7 -f iV, (4) 

and N = n,p, 


accounting for all possible charge combinations. The im¬ 
plementation of the reactions is the same as described 
in Refs. [1, E3- III® other hand, the meson-bbaryon 
processes ® are incorporated here for the first time in a 
transport approach. We describe the relevant cross sec¬ 
tions in section UlI Cl 

4) Photon production in the decays of mesons 
(tt®, 77 , 77 ', w, (j), oi) and the A-resonance, where the parent 
hadrons are produced in baryon-baryon [BB), meson- 
baryon (mB) or meson-meson (mm) collisions in the 
course of the heavy-ion collision. We consider the contri- 
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buttons from the photon decays of the following mesons: 


ttO 

-)> 

7 + 7, 

V 


7 + 7, 

v' 

-)> 

P + 7, 

LO 


7 r° + 7 , 

(t> 

-)> 

7 + 7, 

ai 


TT -1- 7 . 

A 


7 +V, 


The decay probability is calculated according to the cor¬ 
responding branching ratios taken from the latest compi¬ 
lation by the Particle Data Group [1^ , updating slightly 
the value s ap plied in earlier HSD investigations at SPS 
energies [^. The broad resonances - including the 
ai, p, a; mesons - in the initial or final state are treated in 
PHSD in line with their (in-medium) spectral functions 
as implemented and described in detail in Ref. [4l|. 

The photon production from the mesonic decays rep¬ 
resents a ’background’ for the search of the direct pho¬ 
tons. However, this background can only partly be fixed 
by the independent measurements. One usually uses the 
’cocktail’ method to estimate the photon decay spectra 
and their contribution to the elliptic flow V 2 , which re¬ 
lies among others on the mr-scaling assumption for the 
particle spectra. We have assumed throughout that the 
direct photon spectra do not include the contributions 
from the tt, rj, rj' and w decays, because they were sub¬ 
tracted experimentally. 


III. THEORETICAL DEVELOPMENTS 

The extensions compared to our previous publica¬ 
tions on photon production in heavy-ion collisions within 
the PHSD model are described in the following subsec¬ 
tion: an improved implementation of the Bremsstrahlung 
channel m + m —>■ m-l-m-l -7 beyond the SPA is described 
in subsection nil A[ an estimation of the suppression due 
to the Landau-Pomeranchuk-Migdal (LPM) effect will be 
presented in subsection IHI Bl and additional baryonic pro¬ 
cesses will be described in subsection IHI Cl 


A. Bremsstrahlung m-l-m—>^m-|-m -|-7 beyond the 
soft-photon approximation 

In the present work, we improve the description of the 
photon bremsstrahlung in meson-fmeson scattering by 
going beyond the soft-photon approximation [d^. Since 
pions are the dominant meson species in the heavy-ion 
collisions, we concentrate here on the description of the 
bremsstrahlung photon production in pion-l-pion colli¬ 
sions. Indeed, we recall that the pion-|-pion interactions 
are most numerous and provide the dominant source of 
photons from the meson-|-meson bremsstrahlung mech¬ 
anism. This was shown in the previous investiga¬ 


tion 0 explicitly by studying the channel decomposi¬ 
tion of the bremsstrahlung photons. Therefore, we will 
be able to considerably reduce the theoretical uncer¬ 
tainty in modeling the total photon production in the 
hadron bremsstrahlung by improving our modeling of the 
7 r-(- 7 r—)> 7 r-|- 7 r -|-7 process. We will suggest a general¬ 
ization to the other meson species and meson-(-baryon 
interactions in the end of this section. 

In order to calculate the differential cross sections 
for the photon production in the processes of the type 
7r-|-7r —>■ 7r-(-7r-|-7, we use the one-boson exchange 
(OBE) model as originally applied in Ref. [s^ to the 
dilepton bremsstrahlung in pion-|-pion collisions, later on 
in Ref. to the low-energy photon bremsstrahlung in 
pion-fpion and kaon-|-kaon collisions. 

In order to achieve a reliable level of accuracy we model 
the interactions of pions with hadrons using a covari¬ 
ant microscopic effective theory with the interaction La- 
grangian, 


L^nt = gaCrdfj,Tfd'"Tf+gp^-{Tfxdf,n)+gfff,^d'"Tf-d''n, ( 6 ) 

as suggested in Refs. [HjUlj. Within this model the in¬ 
teraction of pions is described by the exchange of scalar, 
vector and tensor resonances: a, p and /2(1270), respec¬ 
tively. Additionally, the form factors are incorporated in 
the vertices in the t- and u-channels to account for the 
composite structure of the mesons and thus to effectively 
suppress the high momentum transfers: 


ha{k'^) 


rnj-ml 


( 7 ) 


where TOq, = Too- or rup or ruf is the mass of the 
exchanged meson and k^ is the momentum transfer 
squared. 

The cross section for tt -|- tt —>■ tt -(- tt scattering is given 
by 


daeijs) ^ lATeiP 

dt 167rs(s — dm^) ’ 

where the matrix element |Mp is calculated by coher¬ 
ently summing up the Born diagrams of the a-, p- and 
/ 2 -meson exchange in t, s and u channels (the n-channel 
diagrams are needed only in case of identical pions), 

IMeip = |M^(a) + M‘(a) + M“(a) 
+M%p)+M\p)+M^{p) 

+M^(/) + M‘(/)+M“(/)|2. (9) 

Let us define the four-momenta of the incoming pions 
as Pa = {Ea,Pa) and pb = (Eb,pb), the momenta of the 
outgoing pions as pi = {Ei,pi) and p 2 = (E 2 ,P 2 ) and 
the four-momentum of the exchanged resonance (tr, p or 
/ 2 ) as k. 

The propagators of the massive and broad scalar and 
vector particles are used to describe the exchange of the 
a and p mesons (see e.g. Ref. [s^). The resonance /2 is a 
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ing relations): 


M‘(a) 

M\p) 

M^[p) 

M\f) 

1 


-glhljt) {2ml-tY 

t — ml + im„T„ 

-gl {s-2mlf 
s — ml + imrjTrj ’ 

_ - 9 phl{t) (s - uf 

t- ml + impTp 

ffp (u - tf 

s — ml + impTp ’ 

t — m?j + imfV f 

- (s - 2 m 2)2 


( 11 ) 


(2m7r^ 



FIG. 1: (Color on-line) Cross section of pion-|-pion elastic 
scattering within two effective models are compared to the 
experimental data from [^, |^: the exchange of two mesonic 
resonances, scalar a and vector p (blue dashed line), and the 
exchange of three resonances a, p and the tensor resonance 
/2(1270) of the particle data booklet (red solid line). The 
green dotted line shows the constant a^i = 10 mb for com¬ 
parison. 


spin-2 particle, for which the full momentum-dependent 
propagator has been derived in [b^. The polarization 
sum is 


PpvaP — padvP + Pfipgva 9)1^901^) 


1 ^ ki/kp 


kvkn 


k,,k. 


9t^p- 






2^1 ^ kpki, ^^1 


■ diyOL 9 

m 
^/3 


/ 


m 


f 


■ 9i^l3 


h h 

( 10 ) 


Following the example of the dilepton production study 
in Ref. [^, we use the same propagator for the /2 res¬ 
onance while additionally accounting for its finite width 
by adding an imaginary part to the self-energy in accor¬ 
dance with the lifetime. 


As the result, the following expressions are obtained 
for the matrix elements of the elastic n + tt scattering 
diagrams (we give here explicitly the t- and s-channel re¬ 
sults, the u-channels can be easily obtained by the cross- 


M%f) 


= g/ 

s — m?^ impTf 

(^^{s - 2miy - {2ml 


tf 


(2m7r^ 



where the Mandelstamm variables are defined as s = 
{Pa + Pbf = {Pl + Pbf, t = {pa - Plf = {Pb - P2f, 
u = {Pa-P2f = {Pb-Plf- 

We point out that the formulae m are compact, 
because the masses of all pions were assumed to be 
equal to m-n and the energy-momentum conservation 
Pa + Pb = Pl + P 2 was assumed. These conditions are 
not satisfied for the off-shell tt + tt ^ tt + ir subprocess, 
which we will encounter in the subsequent calculation of 
the bremsstrahlung photon production 7r-|-7r —>■ Tr-l-Tr-l-y. 
For the actual calculation we will use the off-shell gen¬ 
eralizations M{pa,Pb,Pi,P 2 ) of the formulae (ITT]) , which 
we derive from the Feynman diagrams, but the off-shell 
formulae are too lengthy to be presented here. 

A reduced version of the model with the exchange of 
only two resonances - the scalar a and the vector p me¬ 
son - was used by the authors of Ref. to calculate the 
rate of the photon production from the 7r-|-7r —>■ 7r-|-7r-|-7 
process at low transverse momenta of the photons {pt < 
0.4 GeV). This approximation is suitable at low pr be¬ 
cause the photon rate in this kinematical region is domi¬ 
nated by pion collisions of low center-of-mass energy y/s, 
for which the contribution of the / 2 -exchange is small. 
However, relatively high transverse momenta of photons 
Pt’ = 1 — 2 GeV are of interest for our goal of clarifying 
the ’’puzzling” high elliptic flow of direct photons. Thus 
we need a robust model for the interaction of mesons also 
at y/s > 1 GeV. Therefore, we use the OBE model with 
three mesons as interaction carriers (including the tensor 
particle /2(1270)) in our present calculations. 

Phenomenological coupling constants, masses and 
widths of the three interaction-carriers that enter the La- 
grangian dH]) have to be fixed to the integrated energy- 
dependent cross section of the pion-|-pion elastic scatter¬ 
ing (Tei(s), which is known experimentally. We present 
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COS ( 0 ) 

FIG. 2: (Color on-line) Angular differential cross section for pion-fpion elastic scattering within two effective models: the 
exchange of two mesonic resonances, scalar cr and vector p (blue dashed line), and the exchange of three resonances a, p 
and the tensor particle /2(1270) (red solid line with symbols and gray dashed line). The red solid line with star symbols 
shows the model with the full momentum dependence of the /2 propagator, while the grey dashed line is obtained neglecting 
the momentum dependence of the /2 propagator. The green dashed line shows the constant and isotropic Get = 10 mb for 
comparison. 


in Fig. [T] the integrated cross section of the tt -|- tt elas¬ 
tic scattering in the two versions of the OBE model de¬ 
scribed above: taking into account 2 resonances a, p 
(dashed blue line) and taking into account 3 resonances 
a, p and /2 (solid red line). Fitting the parameters of 
both variants of the OBE model (with two- or three- 
resonance exchange), the data from [^, |6^ can be de¬ 
scribed. The best-fit parameters for the model used here 
are: = 2.0, = 0.525 GeV, F,^ = 0.100 GeV, 

Pp = 6.15, rrip = 0.775 GeV, Fp = 0.15 GeV, gfirif = 8.0, 
ruf = 1.274 GeV, F/ = 0.18 GeV. The values of the 
masses and widths suggest the identification of the p- 
resonance to the p-meson and of the particle /2 to the 
/2(1270) in the particle data book [^ . 

One sees in Fig.[T]that the tensor particle /2 is impor¬ 
tant for the description of the pion interaction at higher 
collision energies ^/s > 1 GeV. Neglecting the contribu¬ 
tion of the /2 leads to an underestimation of the tt -|- tt 
elastic scattering cross section by an order of magnitude 
around y/s = 1.2 — 1.3 GeV. Later data on the n + tt 
interaction at i/s above 1 GeV - extracted in Ref. 
from the measurement of the K +p ^ A-|-7r-|-7r reaction 
- also point to the importance of the tensor interaction 


in the resonance region of the /2(1270). 

With the parameters fixed to the integrated pion 
elastic scattering cross section as described above, we 
can calculate the differential cross section daei{s)/dt, 
which is a necessary ingredient for the calculation of 
the bremsstrahlung photon production within the soft- 
photon approximation for the low photon energy and 
low y/s (see below). The results of our calculations for 
the differential cross section of the process tt + tt ^ tt + tt 
as a function of the scattering angle 0 in the 7r-|-7r center- 
of-mass system is presented in Fig. [2] within two effective 
OBE models: including the exchange of two mesonic res¬ 
onances, scalar a and vector p (blue dashed line), and 
including the exchange of three resonances cr, p and the 
tensor particle /2 (red solid line with symbols and gray 
dashed lines). The red solid line with star symbols shows 
the latter model with the full momentum dependence of 
the /2 propagator, while the grey dashed line is obtained, 
if the momentum dependence of the /2 propagator is ne¬ 
glected. The green dashed line shows the constant and 
isotropic = 10 mb for comparison. 

One can see the pronounced influence of the tensor res¬ 
onance exchange at higher y/s of the collision. One ob- 
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FIG. 3: (Color on-line) Feynman diagrams for photon production in the reaction 7r-|-7r—>'7r + 7r-|-7in the one-boson exchange 
(OBE) model. The time goes from left to right. For identical pions, e.g. the u-channel diagrams have to be added. 


serves that for low ^/s = 0.5 GeV the a peak dominates 
the cross section and the angular dependence is governed 
fully by the combination of scalar and vector terms. On 
the other hand, the influence of the / 2 -resonance is clearly 
seen at -Js = 1 GeV and -/s = 1.5 GeV: the model with 
two resonances (dashed blue line) differs from the results 
of the other OBE variants. The correct momentum de¬ 
pendence of the /2 propagator is important, as one sees 
from the difference of the star red and grey dashed lines. 

Figs. [T] and [2] demonstrate that the inclusion of the 
tensor resonance in the model leads not only to an in¬ 
crease of the cross section for ^/s > 1 GeV (in agreement 
with the data) but also to a considerable ’’flattening” of 
the cos 0 distribution. The cross section at higher ^/s in 
the extended model is more isotropic. 

Using the OBE model (described above) for the covari¬ 
ant interactions of pions, we can calculate the emission of 
photons by the colliding pions by gauge coupling to the 
external hadron lines. The Feynman diagrams for the 
photon production in the process 7r-l-7r—^-Tr-I-Tr-I-T are 
shown in Fig. [S] For identical pions, e.g. 7r+ -b 7r+, the 
u-channel diagrams have to be added, which are obtained 
from the t-channel diagrams by exchanging the outgoing 
pions. The applicability of this method is not limited by 
the low energy of the photon but is restricted only by the 
applicability of the effective model to the description of 
the pion (elastic) interaction. 

Let us again denote the four-momenta of the incoming 
pions by Pa and the momenta of the outgoing pions by 
Pi and p2, while the photon momentum will be denoted 
by q = iqo,q)- The cross section for photon production 


in the process 


7r(po) + 7r(pb) -)> tt { pi ) -b tt { p 2 ) + 7 ( 9 ) (12) 

is given by 

1 


da^ = 


2y^s{s - Ami) 


|M(7)|"di?3, (13) 


where dR^ is the three-particle phase space, which de¬ 
pends on the momenta of the outgoing pions and of the 
photon. 


77 ? - d?‘p 2 d^q 

^ ~ {2Tr)^2Ei {2Rf2E2 {2TTY2qo 

{2RY5^{Pa+Pb-Pl-P2-q)- (14) 

The cross section (|T^ will be integrated over the final 
pion momenta to obtain da/d^q. The (5-function allows to 
perform four integrations analytically and the remaining 
two are done numerically. 

The matrix element M in (jl3|) is a coherent sum of the 
diagrams presented in Fig. [3]- i.e. of the photon attached 
to each pion line tTq, Rb, and '^2 ~ and of contact terms, 
which account for the emission from the vertices and the 
internal lines: 

mi)? = m;{i)m^{i) 

= \M^+Mm ■ ( 15 ) 

The complex matrix elements for the photon emission 
from each of the pion lines M/ are calculated as sums of 





the three meson exchanges (cr, p, / 2 ). For instance: 


with 


Mf = eJ^[M^i{pa,Pb,Pi + q,P2) 

+Mli{pa,Pb,Pi + q,P2) 

+ M^l{Pa,Pb,Pl + q,P2)] , 


'^a,b ~ 


7M 

•^ 1,2 


Qa^b' 

= Qi,2 


( 2 pa^b - qY 
2 pa,b ■ q 
(2pi,2 - qY 
2pi,2 • q 


(16) 


(17) 


where Qi are the charges of the pions in terms of the elec¬ 
tron charge e. The matrix elements for the pion elastic 
subprocess Mei{pa,Pb,Pi + q,P 2 ) are the off-shell gener¬ 
alizations of the formulae CB- „ 

The contact term is taken from Ref. [i^, eq. (14), 
where it was derived by demanding the gauge invariance 
of the resulting cross section. Indeed, the gauge invari¬ 
ance of the result often has to be restored in calcu¬ 
lations within effective models. In the present work, we 
have used the contact terms in order to cancel the gauge- 
dependent parts in the matrix element as in Ref. |65l |. 
Alternatively, one can take into account additional di¬ 
agrams with the emission of photons from the internal 
lines (see Refs [lH) but this method does not always 
eliminate the need for contact terms (see Ref. 0). We 
have verified that q^M^{'^) = 0. Comparing our results 
to calculations with a different gauge-fixing method will 
allow to quantify the uncertainty of the effective model 
applied (work in progress). 

Within the soft-photon approximation (SPA) one 
assumes that the strong interaction vertex is on-shell. 
The soft photon approximation is based on the first-order 
expansion in the Low theorem and is valid at low 
photon energy and low y/Smm of the meson-|-meson col¬ 
lision, as has been studied in detail for the production 
of dileptons in Ref. [s^. In this case the strong inter¬ 
action part and the electromagnetic part can be sepa¬ 
rated, i.e. the soft-photon cross section for the reaction 
mi -I- m 2 —>■ mi -I- m 2 -I- 7 can be written as 


90 


da^{s) 

cPq 


47r^ 


— A(s,r; 


\e-J{q,t)\ 

l)h 


2 d<Jel{s) 


dt 


dt, (18) 


where is the fine structure constant, t is the mo¬ 
mentum transfer squared in the tt -|- tt —>■ n + n sub¬ 
process, and e is the photon polarization. is the 
electro-magnetic current 


.r = -Qa 


■ — Qb 


Pb 


Qi 


Pi 


Q 2 - 


P 2 


{pa ■ q) {pb ■ q) {pi ■ q) {p2 ■ q)' 


The polarization sum 


pol X 


(19) 



FIG. 4: (Color on-line) Cross section for the production of a 
photon with energy go = 0.005 GeV in the process tt -|- tt— > 
TT -f TT -|- 7 within the following models: the exact OBE cross 
section within the effective model taking into account scalar, 
vector and tensor interactions via the exchange of cr, p and 
/2(1270)-mesons gives (red line with star symbols), the soft 
photon approximation to this model (blue dotted line); the 
OBE result within the model taking into account only the 
scalar and vector interactions via the exchange of cr and p 
mesons (blue dashed line), and the soft photon approximation 
to this model (cyan dash-dot-dotted line). 


depends on the photon momentum q, the charges of the 
pions Qi as well as on the invariant kinematic variables, 
including t. For the case of equal-mass particle scattering 
{ma = mb = mi = m 2 = m.^.), one obtains [6l|: 


|e-Jp — — {-iQl +QI + QI + Q2) 
% 


— 2 {QaQb + QlQ 2 ) 
+ 2 {QaQl+QbQ 2 ) 
+ 2 {QaQ 2 + QbQ3) 


s - 2ml + v^s - Ami 


i/s{s - Ami) \ a/s “ ~ J 

2 m‘i — t , f\J—t + 4m2 -j- 


: In 


y/t{t - Ami) \\/-t-b4m2 - 
s — 2mJ -I- t 
y^{s + t){s + t - Ami) 


X In 


Vs + t + ^s + t- Ami 
y/s + t — ^ys + t — Ami 


.( 20 ) 


In (fT51) . d(Jei{s)/dt is the on-shell differential elastic 7r-|-7r 
cross section, which is a function of the invariant energy 
s and the pion scattering angle via t. 

The expression m is considerably simpler in com¬ 
parison to the “exact” OBE formula (TT^ because of the 
factorization of the diagrams from Fig.[3]into the electro¬ 
magnetic part and the elastic 7r-|-7r —>■ 7r-|-7r sub-process, 
for the cross section of which the g-dependence is omit¬ 
ted. This corresponds to neglecting the off-shellness of 
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the pion, which emits the photon, e.g. for the pion a: 

Pa-q^Pa- (21) 

Consequently, the sub-process invariant energy S 2 is also 
approximated by the total invariant energy of the process 

TT + TT^TT + TT + 'y: 

S2 = {Pa+Pb- qf ^ {Pa+Pbf = S, (22) 


and the limits of integration over t are also taken as for 
the on-shell case, i.e. from —A(s, m^, to^)/s to 0, while 
the actual integration over the full 3-particle phase space 
in the exact treatment (TT^ involves different limits for t. 

In Fig. |4] we show the resulting cross sections for the 
photon production in the process 7r-|-7r —>■ Tr-l-Tr-l-y 
within the following models: the “exact” OBE taking 
into account scalar, vector and tensor interactions via 
the exchange of cr, p and /2(1270)-mesons gives the red 
line with star symbols, the soft photon approximation 
(fT^ to this model is shown by the blue dotted line; the 
OBE result within the model taking into account only the 
scalar and vector interactions via the exchange of a and 
p mesons is shown by the blue dashed line, and the soft 
photon approximation to this model is presented by the 
cyan dash-dot-dotted line. The photon energy is fixed 
to go = 5 MeV. For the very low energy of the photon 
of go = 5 MeV the SPA agrees with the “exact” cross 
section very well in the region of ^/s < 0.9 GeV (see 
Fig.lH). However, the discrepancy to the OBE result is 
increasing rapidly with growing -y/s; the calculations for 
the higher photon energy of go = 0.5 GeV show an even 
larger discrepancy between the SPA and the exact OBE 
result. 

Using the cross section for the 7r-|-7r —Tr-l-Tr-l-y 
reaction according to (d as a function of the pho¬ 
ton energy go and the collision energy -y/s, we calcu¬ 
late the yield dN/cPq and the invariant rate qodR/(Pq of 
bremsstrahlung photon production from an equilibrated 
pion gas. Within kinetic theory, the rate of photon pro¬ 
duction in the collisions of particles a and 6 in a ther- 
malized medium (number of photons produced per unit 
space-time volume d'^x) is an integral over the three- 
momenta of the incoming particles: 


qo 


dN 

dx'^d^q 


= g ds 


f ^^Pb (Ea+Eh)/T 

(2^)3 J (27r)3 
da"^ 


>^VreiqQ^^5{s - (Pa+Pbf), (23) 


where T is the temperature, Vrei is the relative velocity 
given by 


'^rel 


{Pa ■ Pb)^ - mlml 

EaEb 


(24) 


and g = {2sa + l)(2st, -I-1) is the spin degeneracy factor. 
Integrating the expression (j23l) over the particle momenta 



% [GeV/c] 


FIG. 5: (Color on-line) Invariant rate of the bremsstrahlung 
photon production from an equilibrated pion gas at a temper¬ 
ature of T = 200 MeV and pion chemical potential pn ~ 0 as 
calculated in the OBE model with three resonance exchange 
within the soft-photon approximation (red dashed line). The 
black solid line ”Haglin:2004” from Ref. is shown for com¬ 
parison. 


one obtains [6l|: 


dN 


qo 


_^ T^g 

d'^xd^q IGtt^ 


■ da^ 

dz --iFi(z)go^, 


(25) 

where Zmin = {'ma + mi,)/T, z = y/sjT, and Ki(z) is the 
modified Bessel function. 

The expression (l23)l can be generalized to account 
for quantum effects such as Bose enhancement or Pauli 
blocking (depending on the particle type) by integrat¬ 
ing additionally over the momenta of the final particles 
and changing the Boltzmann distributions to the Fermi 
or Bose distribution functions fi{T): 


dN ^ f r d^Pa r d^pb f d?pi r d?pi 

dx^d^q J J {2Tr)^ J {2 tt)^ J (27r)3 

fa{T)fb{T){l-MT)){l-f2{T)) 

d(7^ 

XVreiqo-J^d{s -{pa+Pb)'^)- (26) 


In the current section we calculate the thermal rates ac- 
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cording to formula (l23l) . However, within the PHSD 
transport approach for the heavy-ion collisions in sec¬ 
tion |l3 the effects of the quantum statistics will be taken 
into account. 

In Fig. [ 5 ] the rates are presented for a temperature T = 
200 MeV and pion chemical potential p,,r = 0 for the OBE 
model with three resonance exchanges adopting the soft- 
photon approximation (red dashed line). We confirm the 
results from Ref. 0 (black solid line) calculated within 
the same assumptions (SPA, three resonances) but with 
a slightly different parameter set of the Lagrangian. It 
is, however, questionable that the SPA is applicable at 
high photon energies. 

We note that the accuracy of the SPA approximation 
can be significantly improved and the region of its appli¬ 
cability can be extended by slightly modifying the for¬ 
mula m - i-e. by evaluating the on-shell elastic cross 
section at the invariant energy S 2 of the sub-process. The 
latter is kinematically fixed to 

S 2 = s - qo\/s ^ s, (27) 

Thus the modified SPA formula is 


da'^{s) 


a. 


47r 


J k ••/(<?, (28) 




In the following, we will denote the approximation (1281) 
as “improved SPA” and will show below that it provides 
a good description of the exact photon production rates. 

But first we describe a simple model for the pho¬ 
ton production in meson-|-meson collisions by the 
bremsstrahlung mechanism, i.e. the approximation of a 
constant meson + meson elastic cross section. We have 
seen above that the elastic-scattering cross section of pi- 
ons is approaching about cr = 10 mb at high y/s and 
becomes increasingly isotropic. As a very simple esti¬ 
mate for the interaction of two mesons one can use a 
constant isotropic elastic-scattering cross section of, e.g., 
cr = 10 mb at all y/s. In case of the isotropic cross section 
a el, it can be taken out of the integral (fT^ . The integra¬ 
tion of the electromagnetic current over the photon angle 
in m can be done in a straightforward way in case of 
small momentum transfer since 

/” 

J tmin ^ 


This leads to a useful approximation to the cross sec¬ 
tion for the photon bremsstrahlung in meson-f meson col¬ 
lisions 




qo 




cr(s) 


a a(s) 

4^ IT’ 

S-(Mi+M2)2 
- ^172 - 


(30) 

(31) 


where My and M 2 are the masses of the colliding mesons; 
here the first meson is assumed to have unit charge and 


the second meson to be charge-neutral, i.e. the approxi¬ 
mation does not take into account the full dependence on 
the charges of the pions as was done in (l20ll or (ITTll . This 
“constant-cross-section” approximation is useful in par¬ 
ticular for an estimate of the photon bremsstrahlung in 
elastic collisions of mesons, for which the scattering cross 
sections are not known experimentally. We recall that the 
formulae (15(11) and m have been used to model photon 
bremsstrahlung in the collisions of various meson species 
in the previous transport calculations in Refs. [BM. l4ll|. 

Finally, we proceed to calculate the photon production 
rates beyond the soft photon approximation, using the 
cross section for the tt + tt —>■ tt + tt + j reaction calculated 
according to the exact OBE expression (USD. We present 
the calculated invariant rate qodR/dq^ of bremsstrahlung 
photons produced from an equilibrated pion gas at T = 
150 MeV and = 40 MeV in Fig. |6l The results of the 
following models are compared: 

• model 1 (red solid line): exact rates within the 
one-boson exchange model (OBE) beyond the soft- 
photon approximation - i.e. using the formula (1131) 
for the photon production cross section q^ da^/d?q\ 

• model 2 (blue dotted line): result within the soft 
photon approximation - i.e. using the formula (1181) 
- while using the elastic tt -|- tt cross section calcu¬ 
lated within the OBE model as given by the equa¬ 
tions (l5t-(fTT]): 

• model 3 (black short-dashed line): results of the 
improved soft photon approximation - i.e. using 
the formula ((55D instead of (HHD - and the same 
pion elastic scattering cross section as in the model 
2 ; 

• model 4 (green dashed line): soft photon approxi¬ 
mation using a constant isotropic elastic cross sec¬ 
tion of Gel = 10 mb and assuming for the pion 
charges Qa = Qi = 1, Qb = <52 = 0 - i.e. using the 
formula (|5(ID . For this case the elastic cross section 
does not depend on y/s and therefore there is no 
difference between the SPA and improved SPA. 

The rate of bremsstrahlung photons at low transverse 
momenta pT < 0.4 GeV has been calculated before in 
Ref. [15 within the one-boson exchange model with the 
exchange of two resonances for the same system. This 
previous result is shown for comparison by the cyan 
dashed line and is confirmed by our present calculations. 
The agreement is expected, since our calculations differ 
only in the inclusion of the / 2 -meson exchange, which 
is important for large y/s and does not play an impor¬ 
tant role for the production of low transverse momentum 
photons, which is dominated by low y/s of the tt + tt 
collisions. 

On the other hand, the SPA (model 2) deviates from 
the exact OBE result (model I) even at low go: because 
the former directly follows the y/s structure of the elastic 
TTTT cross section. Since the formula (TTSl) does not account 









11 



q„ [GeV/c] 

FIG. 6: (Color on-line) Invariant rate of bremsstrahlung pho¬ 
tons produced from an equilibrated pion gas at T = 150 MeV 
and = 40 MeV versus the photon energy qo- The in¬ 
let shows the same quantity for the range of photon energies 
qo = 0.1 — 0.4 GeV. The calculations have been performed 
within the following models: (1) OBE model beyond the soft- 
photon approximation (red solid line with star symbols); (2) 
OBE model within the soft photon approximation (blue dot¬ 
ted line); (3) OBE model within the improved soft photon 
approximation (black short-dashed line) - the invariant en¬ 
ergy S 2 of the on-shell tt -|- rr elastic process is not equal to the 
total invariant energy of the process s: S 2 = s — qo^/s\ (4) the 
soft photon approximation with the constant isotropic elastic 
cross section of (Jei = 10 using the formula (I30II (green dashed 
line). The cyan solid line ”Liu and Rapp:2007” from Ref. 
is shown for comparison. 


for the off-shellness of the emitting pion, it overweights 
the high-y^ part of the elastic cross section, in line with 
the findings of Refs. [H, [6l|. We note that the OBE 
model presented here is constrained by the pion scat¬ 
tering data only up to = 1.4 GeV and generally 

cannot be extended to large ^/s. Thus the SPA scenario 
’’model 2” is not reliable for large qq (approximately for 
go > 0.8 GeV). This is not the case for the improved SPA 
(model 3). 

One can see in Fig.[6]that the improved SPA (|28|) gives 
a very good approximation to the exact result at higher 
photon energies of up to go ~ 2 GeV. This is because the 


yGJ of the sup-process does not reach such high values 
as ^/s, and the OBE model for the elastic cross section 
is sufficiently realistic in this region of y/sj• 

In comparison, the constant cross-section approxima¬ 
tion (based on formula (1301) ) overestimates the exact rates 
for go > 1 GeV and underestimates for go < 0.4 GeV. 
This model approximately corresponds to the procedure 
used previously in our transport calculations for the es¬ 
timation of the photon bremsstrahlung in meson-fmeson 
collisions in Refs. iami. In the present work, we will 
use the exact OBE cross section da'^/d^q. Therefore, we 
will now find a lower yield of bremsstrahlung photons for 
go > 1 GeV in the transport simulations of heavy-ion 
collisions, see section m below. 

Another reasons for the good agreement between the 
“improved” SPA rates with the exact ones is that the 
dominant contribution to the rates comes from the low 
collision energies ■\/s, while the deviation between the ex¬ 
act cross section of the process 7r-|-7r — ^ 7r-|-7r-|-7 from that 
calculated within the improved SPA is most pronounced 
at high y/s and high go. Such high collision energies y/s 
are suppressed in a thermal medium exponentially by the 
Boltzman factor describing the occupation of the pion en¬ 
ergies at fixed temperature. It is not clear, whether this 
high accuracy of the “improved SPA” holds also out-of- 
equilibrium. We show in Fig. [3 explicitly the distribu¬ 
tions in the number of meson-l-meson, meson-l-baryon 
and parton-bparton collisions versus their invariant en¬ 
ergy y/s as calculated within the PHSD for the case of a 
Au + Au collision with the energy = 200 GeV at 

impact parameter b = 7 fm, which is the centrality qual¬ 
itatively similar to minimum bias. The distribution of 
meson-bmeson, meson-bbaryon and parton-bparton colli¬ 
sion invariant energies dN/dy/s in an actual heavy-ion 
reaction shows a clear dominance of the low--y/s comover 
collisions. Therefore we expect that the deviations of the 
cross section tt + tt ^ Tr-bTr-byin the improved SPA 
from the exact ones at high y/s will not influence much 
the yield of the photons that result from the integration 
over all pion collisions: 

In summarizing the results of the current section, 
we have improved the implementation of the photon 
bremsstrahlung in the process 7r-b7r—^-Tr-bTr-by within 
the PHSD transport approach, using now the exact OBE 
cross sections beyond the soft-photon approximation. 
The bremsstrahlung photon production in collisions of 
other meson types is treated only approximately, i.e. in 
analogy to the tt + ir collisions by means of mass-scaled 
cross sections. 

We note, that another important source of the pho¬ 
tons is the bremsstrahlung in meson+baryon collisions 
(cf. Ref. 0 )- As we have shown above, the SPA 
gives a good approximation to the exact rates, if we use 
the correct invariant energy in the hadronic sub-process 
S 2 = s — qqy/s and a realistic model for the differential 
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FIG. 7: (Color on-line) The distribution in the invariant collision energy ^/s for the elastic scattering of mesons on mesons (left 
panel), mesons on baryons (middle panel) and for the elastic scattering of partons (right panel) in the course of a Au + Au 
collision at ^snn = 200 GeV for & = 7 fm as calculated within the PHSD approach. Here V = (p, a;,<))); K denotes all the 
strange mesons K = [K, K, K*, K*) and B stands for the baryons B = {p,n, A,...). 


cross section of the subprocess, i.e. for the elastic scat¬ 
tering of mesons on baryons. The cross sections for the 
meson-|-baryon elastic scatterings (implemented within 
the PHSD transport approach) have been previously ad¬ 
justed to the data differentially in energy and angular 
distribution. Thus we evaluate the photon production in 
the processes m + B^m + B + j in the PHSD by using 
realistic elastic scattering cross sections taken at the cor¬ 
rect invariant energy in the scope of the improved 
SPA. 


B. The LPM effect 

The radiation of photons by charged particles is modi¬ 
fied in the medium compared to the vacuum. One of such 
medium effects is caused by the absence of well-defined 
incoming and outgoing asymptotic states due to the mul¬ 
tiple scattering of particles in a strongly interacting envi¬ 
ronment. If the subsequent scatterings occur within the 
time necessary for photon radiation Tj ~ 1 /qo, then the 
amplitudes for the emission of photons before and after 
the charged particle scattering have to be summed co¬ 
herently. The effect of this destructive interference on 
the photon spectrum by electrons transversing a dense 
medium was first studied by Landau and Pomeranchuk 
in Ref. [t^ and Migdal in Ref. [^. Accordingly, the 
Landau-Pomeranchuk-Migdal (LPM) effect modifies the 
spectrum of photons produced in the medium in com¬ 


parison to the incoherent sum of emissions in quasi-free 
scatterings, leading especially to a suppression of the low 
energy photons because the formation time of the photon 
Tj is proportional to the inverse photon energy 1/qo- In 
particular, the LPM effect regularizes the 1/qo divergence 
of the quasi-free bremsstrahlung spectra. The LPM sup¬ 
pression and the induced thermal mass of the medium 
quanta (the dielectric effect) together ensure that the 
photon spectrum is finite in the limit go 0. 

The importance of the LPM effect for the case of 
dilepton and photon production from QCD systems was 
shown in Refs. The magnitude of the LPM sup¬ 

pression is governed by the average time between the col¬ 
lisions T, which in turn is given by the inverse scattering 
length a or by the inverse average spectral width of the 
particles P: 


The LPM suppression is more pronounced in case of 
small T, i.e. for high reaction rates. Thus we expect 
it to be important for the emission of photons from the 
strongly-interacting quark-gluon plasma (sQGP) as cre¬ 
ated in the early phase of the heavy-ion collision. In¬ 
deed, it was shown in Refs. [H, 111 in the scope of the 
DQPM that the average collision time of partons is as 
short as r Ri 2 — 3 fm/c for temperatures in the range 
T= 1—2 Tc, where Tc ~ 158 MeV is the deconfinement 
transition temperature. In comparison, the average time 
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% [GeV/c] 

FIG. 8: (Color on-line) Invariant rate of photons produced 
from the strongly-interacting quark-gluon plasma consisting 
of massive broad quasi-particle quarks and gluons (red solid 
line). The leading-order pQCD rate (blue dashed line) from 
Ref. is shown for comparison. 


between pion collisions in a thermalized pion gas at tem¬ 
peratures T <Tc is above 10 fm/c, see Ref. (74 |. 

The photon production in the sQGP proceeds by the 
interactions of dressed quarks and gluons through quark 
annihilation and gluon Compton scattering processes 
g -b g —?> g -b 7 , q{q) + g —t q{q) -b 7 as well as through 
the quark bremsstrahlung and off-shell parton decays. In 
the strongly interacting QGP (and in PHSD), the gluon 
and quark propagators differ significantly from the non¬ 
interacting propagators such that bare production am¬ 
plitudes can no longer be used [13, 113 ■ The off-shell 
quarks and gluons have finite masses and widths, which 
parametrize a resummed interaction of the QGP con¬ 
stituents. The quark off-shellness leads to higher twist 
corrections (^ w^/s, m^/t, m i/u ) to photon and dilep¬ 
ton production cross sections [m . [t^ . These corrections 
are small in hard hadron scattering at high center-of- 
mass energy ^/s > 10 GeV, but become substantial for 
photon production in the sQGP, where the characteris¬ 
tic y/s of parton collisions is of the order of a GeV (cf. 
r.h.s. of Fig. [71). Using the cross sections for photon 
radiation by dressed quarks and gluons in the processes 
qq —>■ 57 and qg —>■ 57 as derived in Ref. , we have cal¬ 
culated the rates of photons produced in a thermalized 
strongly interacting QGP. Fig. [ 8 ] presents the invariant 
rate of photons produced from the strongly-interacting 
quark-gluon plasma at the temperature T = 200 MeV 
(red solid line). The leading-order Log-resummed per¬ 
turbative QCD rate (blue solid line) from Arnold, Moore 
and Yaffe (AMY rate) is taken from Ref. [t^ and shown 


for comparison. One observes a qualitative agreement 
between the results of both models although the degrees- 
of-freedom and their couplings are different. We men¬ 
tion that photon rates recently calculated at the NLO in 
perturbative QGD [sol - l^ also approximately are in line 
with those presented in Fig. El 

Let us now quantify the magnitude of the LPM effect 
on the spectrum of photons radiated from the QGP as 
calculated within the PHSD. The coherent photon pro¬ 
duction rate - taking into account the LPM effect - dif¬ 
fers from the non-coherent cross section by a suppres¬ 
sion factor, which generally depends on the photon en¬ 
ergy, temperature and the interaction strength of the con¬ 
stituents. The coherent photon emission rate was derived 
in Ref. [t^ for an elastically interacting pion gas in the 
soft photon approximation for the photon radiation am¬ 
plitudes. The authors of Ref. [t^ used the same method 
for the calculation of the photon emission over the whole 
trajectory of the charged particle as was adopted in the 
original work by Migdal in Ref. [t^ . After averaging over 
the times between collisions assuming an exponential 
distribution 


dW 


-fa 


(34) 


the coherent photon emission rate was found to be 


dR ^ / .2 ( 1 -cos^e) \ 

dq^ (27r)2 \ -b go(l — w cos 0 )^ / ’ 


(35) 


where the brackets < . > stand for an average over the 
velocities after the scattering (w, cos0), while N is the 
number of scatterings and uem ~ 1/137. A realistic 
parametrization of the data was used for the pion elas¬ 
tic scattering cross section (cf. section IIII Al) , but the 
scattering was assumed isotropic. The incoherent rate is 
obtained from (I35|) in the limit a = 0. 

An analytical form of the coherence factor was ob¬ 
tained in Ref. [t^ in the model of hard scattering centers, 
using a quantum mechanical approach to coherently sum 
the photon amplitudes from all the scatterings. In the 
thermal medium the spatial distribution of the scattering 
centers is assumed random. Consequently, the function 
(1341) naturally arises in this model for the distribution 
of times between collisions by a direct calculation of the 
two-particle correlation function. The quenching factor 
in the dipole limit (g = 0 ) was found as 


(G(gor))^ 


( {gpr? Y 
\l + {Qorr) 


(36) 


Although formula (l36l) was obtained in a simple model, 
it is useful because it correctly captures the dependence 
of the LPM suppression on the average strength of the 
interaction given solely by the mean-free-time between 
collisions r in the assumption of isotropic collisions. 

The perturbative interaction of quarks and gluons is 
dominated by small scattering angles due to the massless 
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particle exchange in the t-channel diagrams. In this case 
the coherence factor for the quark system in the limit of 
small scattering angles was obtained in Ref. [t^. How¬ 
ever, up to now the LPM effect in case of a strongly inter¬ 
acting QGP with dressed broad quasi-particles has not 
been evaluated. The elastic scattering of dressed quarks 
in the PHSD is not dominated by the t —?> 0 pole as in the 
perturbative case since the gluon mass (of order 1 GeV) 
acts as a regulator in the amplitude. Accordingly, the 
angular distribution for quark-quark scattering is closer 
to an isotropic distribution for low or moderate y/s in 
accordance with the model assumptions of Ref. [t^ such 
that the expression (l36l) should apply as an estimate of 
the LPM suppression for the photon emission within the 
PHSD. 


In Fig. Owe show the photon emission rate in a QGP at 
the temperature T=190 MeV as calculated in the PHSD 
as an incoherent sum of the photon emission in quark and 
gluon scatterings (red solid line). The blue dashed line 
gives the same rate with the quenching factor (1361) ap¬ 
plied using t(T) = I/r(T) Ri 3.3 fm/c from the DQPM 
(for T= 190 MeV). We observe that the suppression in 
comparison to the incoherent rate is visible only for pho¬ 
ton energies qo < 0.4 GeV. For an estimate of the up¬ 
per limit on the LPM suppression we employ the relax¬ 
ation time approximation for the ratio of the shear vis¬ 
cosity over entropy density y/s which gives rj/s r; 0.14 
at T=190 MeV in DQPM The lowest bound 

as conjectured within the AdS/GFT correspondence is 
rj/s = l/(47r) Ri 0.08. In the relaxation time approxima¬ 
tion this corresponds to a lower value of r r; 1.9 fm/c. 
The coherent photon rate in this case is given by the ma¬ 
genta dotted line and shows a peak in the photon rate 
for go ~ 0.2 GeV. 

In order to further clarify the strength of the LPM 
suppression of the photon emission in the sQGP, we use 
the knowledge of the electric conductivity cro{T) of the 
sQGP from the DQPM which is roughly in line with 
more recent results from lattice QGD. For details of the 
latter analysis on the electric and magnetic response of 
the sQGP we refer the reader to Ref. [s^. We recall 
that the photon emission rate from a thermal medium is 
controlled by uo via the relation. 



FIG. 9: (Color on-line) Incoherent invariant photon produc¬ 
tion rate from the strongly-interacting quark-gluon plasma 
consisting of massive broad quasi-particle quarks and gluons 
(red solid line) scaled by Stt®/ (ST^) in order to match the elec¬ 
tric conductivity for go 0 (cf. Eq. (1371) 1. The blue dashed 
line and the magenta dotted line show the coherent rates with 
the two assumptions for the average time between the colli¬ 
sions r - from the DQPM model and from the AdS/CFT 
correspondence. 


Taking into account some uncertainty in the determi¬ 
nation of T and the expression (1361) . we conclude from 
Fig. [S] and analogous calculations at different temper¬ 
atures that the LPM effect influences the photon pro¬ 
duction from the QGP at photon energies below go ~ 
0.4 GeV, but is negligible for higher photon energies. We 
note in passing that the suppression of the photon spec¬ 
trum in the hadronic phase is much smaller due to the 
lower interaction rate, i.e. longer interaction time r. 


— = — r f IS?) 

T 3T^ ijo^Su (PqJ ’ ' Additional meson-nucleon processes 


where T is the temperature of the system, go is the pho¬ 
ton energy and g is the photon momentum. Using the 
number for ctq/T from the PHSD at the temperature of 
T = 190 MeV from Ref. [^, we obtain a limiting value 
for the scaled photon emission rate of 0.04 for go —)> 0 
according to formula (1571) (green short dashed line in 
Fig. E]). The blue dashed line in Fig. [5]- the estimate 
of the rate based on formula (l36l) and the DQPM aver¬ 
age spectral width of the quarks/antiquarks - indeed ap¬ 
proaches the limiting value of 0.04 as given by the kinetic 
calculations of the electric conductivity in Refs. 


In the present work we incorporate into the PHSD ap¬ 
proach additionally the 2 —>■ 2 processes V + N ^ A^ + 7, 
where V stands for a vector meson while N denotes a pro¬ 
ton or neutron. These processes are the baryonic counter¬ 
parts to the mesonic 2 —)> 2 reactions tt + pj-K ^ ^ + t: jp. 
We here consider the interaction of nucleons with the 
mesons V = p, (j), uj, taking into account the various 
possible charge combinations, e.g. p'^ + p —>■ j + p, 
p~ + p ^ j + n, p~^ + n ^ j + p, etc. Additionally, 
we take into account the photon production in the decay 
of the A-resonance, A —V -|- 7. 
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In order to evaluate the probabilities for photon pro¬ 
duction in the collisions of vector mesons with nucle¬ 
ons, we use the inverse processes 7 -I- iV —>■ p + N^ 
7 -I--/V -^(l) + N,^ + N^uj + N (controlled by data) and 
employ detailed balance to obtain the differential cross 
sections for the processes p+N —>• 7 -I--/V, (j)+N —>■ 7 -I- 
ca -j- N —y 7 “f A^- 

We recall that the detailed balance formula reads: 

*2 

a{NV^jN) = i:LE^ cr{jN ^ NV), (38) 
9v Pnv 

where g-y = 2 and gv = 5 are the spin degeneracy factors 
of the photon and the vector meson V, p%y is the center- 
of-mass momentum in the A^ -I- 7 system, and p*j^y is the 
center-of-mass momentum in the A^ -I- system. 

The cross sections for the exclusive photo-production 
of p, (j) and oj vector mesons on the nucleon have 
been measured by the Aachen-Berlin-Bonn-Hamburg- 
Heidelberg-Munich (ABBHHM) Collaboration and pub¬ 
lished in Ref. [s^. In the same work also parametriza- 
tions for these cross section have been given that are 
based on the vector-meson-dominance model with a non- 
relativistic Breit-Wigner (BW) spectral function for the 
p-meson. Later, the fits have been updated in Ref. [s^ 
using relativistic BW spectral functions for p, w and (j) 
mesons. 

The total cross sections - fitted in Ref. [s^ to the data 
from Ref. [s^ - are given by 

a{-fN ^ VN) = ^— f dp\MvfpNv^vip), (39) 
Pnj^ j 

where the mass of the vector meson is distributed by the 
spectral function Av{p): 


^ 2_ M^r(p) _ 

TT - Mfy + 


(40) 


with Mi denoting the pole mass of the meson. The matrix 
elements for 7 -|- A^ —>■ R -|- A^ are parametrized as 




0.16 mb GeV^ 


0.08 Pl^N 


2(yi-1.73 GeV)2+p^2^ 
0.004 mb GeV^ 


mb GeVR 


(41) 


The cross sections (|3^ with the parameters (ITO are con¬ 
sistent with the dynamics of vector mesons in the PHSD, 
where also relativistic BW spectral functions for vector 
mesons are used and propagated off-shell. 

For the angular distribution of the p-meson production 
in the process 7 -I- A^ —>■ A^ -f p, we follow the suggestion 
of Ref. [1^ , 


da 

dt 


exp(Rt), 


(42) 


with the photon-energy dependent parameter B (fitted to 
the data): B = 5.7 for qq < 1.8 GeV, B = 5.43 for 1.8 < 


go < 2.5 GeV, B = 6.92 for 2.5 < go < 3.5 GeV, B = 8.1 
for 3.5 <go< 4.5 GeV, B = 7.9 for go > 4.5 GeV. 

The data in Ref. [s^ have shown that the cross sec¬ 
tion is dominated by the t « 0 region in line with 
the physics assumptions of the vector dominance model 
(VDM) where the process 7 + V—^-R-I-Vis described 
by the incident photon coupling to the vector meson of 
helicity ± 1 , which consequently is scattered elastically by 
the nucleon (cf. Refs. [stI - I^ I. 

Let us now briefly describe the modeling of the pho¬ 
ton production in the decays of the A-resonance. The 
A —>■ N'y width depends on the resonance mass Ma, 
which is distributed according to the A spectral func¬ 
tion. Starting from the pioneerin g w ork of Jones and 
Scadron [^, a series of models |9ll - l^ provided the 
mass-dependent electromagnetic decay width of the A 
in relation to the total width of the baryon. We employ 
the model of Ref. in the present calculations where 
the spectral function of the A resonance is assumed to 
be of relativistic Breit-Wigner form. We adopt the ” Mo- 
niz” parametrization for the shape of the A spectral 
function, i.e. the dependence of the width on the mass 
r*°*(MA). 


IV. COMPARISON TO DATA FROM 
HEAVY-ION COLLISIONS 

Before presenting the results on photon spectra from 
heavy-ion collisions we stress again that the PHSD ap¬ 
proach so far has provided a consistent description of the 
bulk properties of heavy-ion collisions - rapidity spec¬ 
tra, transverse mass distributions, azimuthal asymme¬ 
tries of various particle species - from low SPS up to LHC 
energies . Furthermore, dilepton production 

from hadronic and partonic sources has been calculated 
at SPS, RHIC and LHC energies and success¬ 

fully compared to the available dilepton data [9^. Ac¬ 
cordingly, the global dynamics of the bulk matter, its 
collective flow as well as the electromagnetic emissivity 
in heavy-ion collisions appear to be well under control. 
Here we extend our previous photon studies in Refs. @,0] 
at the top RHIC energy also to the top SPS and LHC 
energy incorporating the improved cross sections as de¬ 
scribed above. 

The inclusive photon yield as produced in p + p and 
A + A collisions is divided into “decay photons” and “di¬ 
rect photons”. Decay photons - which constitute the 
major part of the inclusive photons - stem from the pho¬ 
tonic decays of hadrons (mesons and baryons). These 
decays occur mainly at later times and outside of the 
active reaction region and therefore carry limited infor¬ 
mation on the initial high-energy state. Consequently, 
it is attempted to separate the decay photons from the 
inclusive yield (preferably by experimental methods) and 
to study the remaining, direct photons. One usually uses 
the “cocktail” method to estimate the contribution of 
the photon decays to the spectra and to the elliptic flow 
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V 2 , which relies among others on the mT-scaling assump¬ 
tion and on the photon emission only by the finally pro¬ 
duced hadrons with the (momenta) distributions as in 
the final state. Depending on the particular experimen¬ 
tal set-up, different definitions of the decay photons are 
applied: whether only the decays of tt^- and ? 7 -mesons 
are attributed to the decays photons or also the decays 
of the less abundant and shorter-living particles 77 ', w, (p, 
ai and the A-resonance. Indeed, the determination of 
the latter contributions (in particular, oi. A) by experi¬ 
mental methods is questionable, because of the emission 
during the absorption and regeneration in the initial in¬ 
teracting phase. Therefore, the theoretical understand¬ 
ing of the decay photon contributions to the inclusive 
spectrum is important. Especially when analyzing simul¬ 
taneously various measurements at different energies and 
within different experimental settings. One can consider 
theoretically decays photons from the following processes: 

7 r° —>• 7 -I- 7, 77 —>■ 7 + 7 , 77' —^ p -I- 7, 

w —>■ 7r° -I- 7, 0 —>■ 77 -I- 7, oi — >■ TT -I- 7, A —>■ 7 -I- A, 

where the parent hadrons are produced in baryon-baryon 
{BB), meson-baryon (mB) or meson-meson (mm) colli¬ 
sions in the course of the heavy-ion collision. 

The following contributions to the direct photons have 
been identified so far: 

• The photons at large transverse momentum pT, so 

called “prompt” or “pQCD” photons, are produced 
in the initial hard N + N collisions and the from 
the jet fragmentation, which are well described the 
perturbative QCD (pQCD) calculations. The lat¬ 
ter, however, might be modified in A + A contrary 
to p+p due to the modification of the parton distri¬ 
butions (initial state effect) the parton energy loss 
in the medium (final state effect). In the A + A 
collisions at large px, there may also be contri¬ 
butions from the induced jet-7-conversion in the 
QGP and the jet-medium photons from the scat¬ 
tering of hard partons with thermalized partons 
qhard + q(g)QGP 1 + qig), but these are relatively 

small. The prompt photons are well modeled by 
using the perturbative calculations and adjusting 
them to the high-p^ region of the observed direct 
photon yield. 

• After the deduction of the prompt photons from 
the direct photon spectra, there is a significant re¬ 
maining photon radiation at pr < 3 GeV observed, 
which is dubbed as “thermal” photons. The low-pr 
photons are emitted by the various partonic and 
hadronic sources as listed below: 

1. Photons are radiated by quarks in the inter¬ 
action with other quarks and gluons, such as 
in the reactions 

q + q^ g + i, q/q + g^q/q + i- 


In addition, photon production in the 
bremsstrahlung reactions q+q/g ^ q+q/g +1 
are possible. 

2. All colliding hadronic charges (meson, 
baryons) can also radiate photons by the 
bremsstrahlung processes: 

m + m^m + m + 'y m + B^m + B + "f. 

3. Additionally, the photons can be produced in 
binary meson-|-meson and meson-|-baryon col¬ 
lisions. We consider within the PHSD the di¬ 
rect photon production in the following 2 —>■ 2 
scattering processes: 

7r-|-7r—^-p-l-y, tt + p ^ n + 

V + N ^j + N, 

where V = p, cp, w, and N = n,p, 

accounting for all possible charge combina¬ 
tions. 


A. Direct photon spectra from SPS to LHC 

We start with the system Pb-|-Pb at ^JsNN = 17.3 
GeV, i.e. at the top SPS energy. Fig. [in] shows the 
comparison of the PHSD calculations to the data of the 
WA98 Gollaboration from for 10% centrality in the 
pseudorapidity interval 2.35 < 77 < 2.95. In addition to 
the sources, which have been incorporated in the origi¬ 
nal HSD study [4l|, the meson-|-baryon bremsstrahlung, 
W —>■ A 7 , A —>■ Nj decay and the QGP channels 
are added. Gompared to the earlier results of Ref. (4ll |. 
the description of the data is further improved and the 
conclusions remain unchanged: the bremsstrahlung con¬ 
tributions are essential for describing the data at low 
Pt- This interpretation is shared by the authors of 
Refs. [^,[ 711193 , who also stressed the importance of the 
meson-l-meson bremsstrahlung in view of the WA98 data 
using hydrodynamical or fireball models. Note that the 
photon contribution from the QGP is practically negligi¬ 
ble for low Pt and reaches at most 25% at px > 0.5 GeV. 

We now step on to the top RHIG energy of y/sNN 
— 200 GeV and update our results on the differential 
photon spectra for the system Au-I-Au. The results for 
the direct photon spectrum as a sum of partonic as well as 
hadronic sources for the photons produced in minimum 
bias Au-|-Au collisions is presented in Fig.lTlIas a function 
of the transverse momentum at mid-rapidity \y\ < 0.5. 
This observable has been calculated within the PHSD 
earlier ig, but has to be revisited in the present work 
since we have incorporated additional photon production 
channels such as the binary baryon-fmeson collisions V -I- 
N —>■ A-I- 7 , A—>■ A 7 and improved the calculation 
of the m + m —>■ 77 ^- 1 - 771 - 1-7 bremsstrahlung channel 
based on the OBE model results which go beyond the 
soft-photon approximation. Indeed, a direct comparison 
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[GeV/c] 


FIG. 10: (Color on-line) Comparison of the PHSD calcu¬ 
lations to the data of the WA98 Collaboration from [^ . 
In comparison to the original HSD study |4i|: (i) the me- 
son-|-baryon bremsstrahlung (blue dash-dotted line), A de¬ 
cays (black short-dashed line) and the photons from QGP 
(green line with round symbols) are added (ii) and the me- 
son-|-meson bremsstrahlung is now calculated beyond the SPA 
(magenta dashed line). The black line with diamond symbols 
labeled as “other” includes: w, r]', (f> an d oi-meson decays, 
binary channels tt -|- pj-K -k jp 7 and A -|- P —>■ A -f 7 . 


of Fig. [TT] to our previous results for the direct photon 
spectra at RHIC (cf. Fig. 1 in Ref. [ 3 ] and Fig. 3 of 
Ref. Q) shows that the data seem to favor the OBE 
model employed here in comparison to the SPA with the 
simple constant isotropic cross section employed earlier. 
In particular the agreement with the data is improved in 
the high-pT range. 

We recall that in our calculations the direct pho¬ 
ton spectrum has the following contributions: photon 
bremsstrahlung in meson-|-meson m + m—+ y 
(blue dash-dotted line) and meson-|-baryon collisions 
m + B —^ m + B + ^ (magenta dashed line); photon 
production in the QGP in the processes q + q ^ + 7 ) 

and q{q) + g ^ q{q) + 7 (green line with yellow circles); 
the reactions tt + p —>■ tt -|- 7 and tt + tt —>■ p -|- 7, (cyan 
line with open symbols); the meson-fbaryon reactions 
P + A —>■ A -b 7 (orange line with filled symbols); decays 
of (j) and oi mesons and of A-baryons (black short-dashed 
line); the photon production in the initial hard collisions 
(”pQCD”) is given by the hard photon yield in p-bp col¬ 
lisions scaled with the number of collisions Ncoii (violet 
dotted line). The measured transverse momentum spec¬ 
trum dN/dpT (given by the filled circles) is reproduced 



Pt [GeV/c] 

FIG. 11: (Color on-line) PHSD results for the spectrum of di¬ 
rect photons produced in minimum bias Au-bAu collisions at 
-^snn = 200 GeV as a function of the transverse momentum 
Pt at mid-rapidity \y\ < 0.5. The data of the PHENIX Col¬ 
laboration are taken from Refs. [^, |^. For the individual 
lines see the legend in the figure. 


well by the sum of all partonic and hadronic sources (red 
solid line). 

As we have previously stressed in Ref. Q, the cen¬ 
trality dependence of the photon spectra carries addi¬ 
tional information, which can be used to disentangle the 
hadronic and partonic sources of the photon produc¬ 
tion. The recent measurements by the PHENIX Col¬ 
laboration O confirm the predictions within the PHSD 
from Ref. 3]- The centrality dependence of the inte¬ 
grated thermal photon yield in PHSD was found to scale 
as A“^^( with the exponent a ~ 1.5, which is in a 
good agreement with the most recent measurement of 
a = 1.48±0.08±0.04 by the PHENIX Collaboration [HI. 
We compare our calculations and the data in Fig. [121 In 
the present work we have added the cross sections of the 
processes F-bA—^-y-bAby exploiting detailed bal¬ 
ance from the cross sections (1551) and the photons from 
the A-Dalitz decays. The contributions of these channels 
to the direct photon yield in the Au -b Au collisions at 
y/SNN = 200 GeV at four different centralities are ex¬ 
plicitly shown in Fig. 1121 where the recent measurements 
by the PHENIX Collaboration (HI are compared to the 
PHSD spectra, too. It can be seen in Fig. [T5] that the 
inclusion of the photon production from the F -b A re¬ 
actions does not enhance the total direct photon yield 
much at all four centralities. We find that partonic chan¬ 
nels and bremsstrahlung photon production “overshine” 
the 2 —>■ 2 reactions, both meson-bmeson and (vector) 
meson-bbaryon channels. 









































18 



Pt [GeV/c] 


Pt [GeV/c] 


FIG. 12: (Color on-line) Contribution of the photon production in the two-to-two p-|-nucleon interaction (orange dashed lines) 
to the total direct photon spectra (red lines) at the top RHIC energy for different centralities. The dominant sources are 
the photons from the QGP and from the hadronic two-to-three bremsstrahlung processes. The PHENIX data are from Refs. 


the phot 

[ 13 , 0 . 


The modification of vector meson spectral functions in 
the nuclear medium is an interesting phenomenon, which 
has been found essential for the understanding of dilepton 
production at SPS and RHIC energies [5lj,[53. Accord¬ 
ingly, there should be also an effect of the in-medium 
modification of vector-meson spectral functions on the 
photon production in the V + N reactions. However, 
since the V + N channel does not dominate the spec¬ 
trum at any photon energy, we conclude that this effect 
cannot be resolved in the total direct photon spectra at 
the present level of the data accurac y a t RHIC; this has 
been explicitly shown earlier in Ref. [4l| for the top SPS 
energy. 

Next, let us investigate the photon production across 
the phase transition in the heavy-ion collision to check 
whether the observed yield of direct photons is produced 
dominantly in some particular region of the energy- 
density or in some particular phase of matter. Fig. 1131 


shows the yield of photons produced at mid-rapidity in 0- 
20 % most central Au+Au collisions at ^/snn = 200 GeV 
as functions of the approximate local “temperature” (i.e. 
the fourth-root of the energy density) from the PHSD. 
The top panel of Fig, [represents the calculations for pho¬ 
tons with the transverse momentum px = 0.5 GeV, while 
the bottom panel corresponds to photons with a trans¬ 
verse momentum px = 1.5 GeV. We observe that the 
early, hot state does not dominate the photon production 
in the QGP contrary to expectations of the static ther¬ 
mal fireball model, where photon production is roughly 
proportional to a power of the temperature (^ T"^). The 
integration over the dynamical evolution of the heavy- 
ion collision leads to roughly the same contribution of 
the different energy density regions since the rate de¬ 
creases but the space-time volume increases. The pho¬ 
ton production in the hadronic phase is dominated by 
the lower energies/temperatures because of the very long 
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FIG. 13: (Color on-line) Top: Yield of photons with the 
transverse momentum pr = 0.5 GeV at mid-rapidity pro¬ 
duced in 0-20 % most central Au Au collisions as func¬ 
tions of the approximate local “temperature” (the fourth-root 
of the energy density) from the PHSD from meson-meson 
bremsstrahlung (dash-dotted lines) and gluon Compton scat¬ 
tering (solid lines). Bottom: Same as in the top panel for 
photons with the transverse momentum pr ~ 1.5 GeV. 


times over which the produced hadrons continue to in¬ 
teract elastically, which is accompanied by the photon 
bremsstrahlung in case of charged hadrons. 

We now increase the invariant energy y/SNN by 
roughly another factor of 14. In Fig. [T3] we show the 
direct photon yield from PHSD in Pb-|-Pb collisions at 
the invariant energy ^/snn = 2.76 TeV for 0-40% cen¬ 
trality. When comparing the theoretical calculations to 
the preliminary data of the ALICE Collaboration from 
Ref. 0 , we find a rather good overall agreement with the 
data within about a factor of 2 in the range of transverse 
momenta pp from 1 to 4 GeV. On the other hand, the 
calculations tend to underestimate the preliminary data 
in the low-pT region 

Furthermore, we provide predictions for the centrality- 
dependence of the direct photon transverse momentum 
distributions at the LHC energy in Fig. [HI since the dif¬ 
ferential investigation in centralities and energies of the 



FIG. 14: (Color on-line) Yield of direct photons at mid¬ 
rapidity in Pb+Pb collisions at the invariant energy y/SNN = 
2.76 TeV for 0-40% centrality within the PHSD in comparison 
to the preliminary data from the ALICE Collaboration Q|. 


heavy-ion collisions will provide crucial information for a 
clarification of the relative importance of the contribut¬ 
ing photon sources. Very recently, the data of the ALICE 
Collaboration for the photon yield at these three centrali¬ 
ties has become available [9^ , we show them as the black 
symbols in Fig. [151 

In conclusion, we have found that from SPS to LHC en¬ 
ergies the radiation from the sQGP constitutes less than 
half of the observed number of direct photons for central 
reactions. The radiation from hadrons and their interac¬ 
tion - which are not measured separately so far - give 
a considerable contribution especially at low transverse 
momentum. The dominant hadronic sources are the me¬ 
son decays, the meson-meson bremsstrahlung and the 
meson-baryon bremsstrahlung. While the former (e.g. 
the decays of w, rj’, (j) and ai mesons) can be subtracted 
from the spectra once the mesonic yields are determined 
independently by experiment, the reactions 7r-|-p —>■ Tr-l-y, 
TT + TT ^ p + V + N ^ N + ^, A—j-fV-l-yas well 
as the meson-meson and meson-baryon bremsstrahlung 
can be separated from the partonic sources only with the 
assistance of theoretical models. 


B. Elliptic flow of direct photons 


The azimuthal momentum distribution of the emitted 
particles is commonly expressed in the form of a Fourier 
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FIG. 15: (Color on-line) Transverse momentum spectra of direct photons at mid-rapidity in Pb-|-Pb collisions at the invariant 
energy ^/snn = 2.76 TeV for the three centrality bins - 0-20%, 20-40% and 40-80% - as predicted within the PHSD. The 
black symbols show the very recently available data from the ALICE Collaboration [^. See the legend for the individual 
contributions. 



FIG. 16: (Color on-line) Elliptic flow V 2 versus transverse 
momentum pr for the direct photons produced in minimum 
bias Au-I-Au collisions at ^/sNN = 200 GeV calculated within 
the PHSD (solid red line); the blue band reflects the uncer¬ 
tainty in the modeling of the cross sections for the individual 
channels. The data of the PHENIX Collaboration are from 
Ref. [H. 

series as 

d^N d^N ( / i T 

' 2.prdpTd, (l+i:2"n(PT)co.[n» - r.„)]j , 

where Vn is the magnitude of the n'th order harmonic 
term relative to the angle of the initial-state spatial plane 
of symmetry and p = {E^p) is the four-momentum 


of the particle under consideration. We here focus on 
the coefficients V 2 and which implies that we have 
to perform event by event calculations in order to catch 
the initial fluctuations in the shape of the interaction 
zone and the event plane '^ep- We calculate the trian¬ 
gular flow U 3 with respect to 'F 3 as U 3 {'I' 3 } = (cos(3[^ — 
'I' 3 ]))/Res('I' 3 ). The event plane angle if 3 and its res¬ 
olution Res(l[' 3 ) are evaluated using the hadron-ha dron 
correlations at larger rapidities as describ ed in Ref. |infll | 
via the two-sub-events method 

We recall that the second flow coefficient V 2 carries 
information on the interaction strength in the system - 
and thus on the state of matter and its properties - at 
the space-time point, from which the measured particles 
are emitted. The elliptic flow V 2 reflects the azimuthal 
asymmetry in the momentum distribution of the pro¬ 
duced particles {px vs Py), which is a consequence of the 
geometrical azimuthal asymmetry of the initial reaction 
region. If the produced system is a weakly-interacting 
gas, then the initial spatial asymmetry is not effectively 
transferred into the final distribution of the momenta. 
On the contrary, if the produced matter has the proper¬ 
ties of a liquid, then the initial geometrical configuration 
is reflected in the final particle momentum distribution. 


More than a decade ago, the WA98 Collaboration 
has measured the elliptic flow V 2 of photons produced 
in Pb + Pb coll isions at the beam energy of Ebeam = 
158 AGeV [Tol, and it was found that the ^ 2 ( 7 ““') 
of the low-transverse-momentum inclusive photons was 
equal to the ^ 2 ( 7 ’^) of pions within the experimental un¬ 
certainties. This observation has lead to the conclusion 
that either (Scenario 1:) the elliptic flow of the direct 
photons was comparable in magnitude to the ^ 2 ( 7 *"''^*) 
and ^ 2 ( 7 '^““^), or (Scenario 2:) the contribution of the 
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t [fm/c] 


FIG. 17: (Color on-line) Rate of photon production in the 
QGP versus time for 0-20% Au-|-Au collisions at y^sjvjv = 200 
GeV from the PHSD. The insert shows the time evolution of 
the elliptic flow of photons from the partonic sources. 



FIG. 18: (Color on-line) Elliptic flow V 2 versus transverse mo¬ 
mentum pt for the charged particles produced in 20-30% cen¬ 
tral Pb-|-Pb collisions at ^sjvjv = 2.76 TeV from the PHSD 
(solid red line ) in comparison to the data from the ALICE 
Collaboration [l03l] : the blue band reflects the statistical un¬ 
certainty of the PHSD calculations. 


direct photons to the inclusive ones was negligible. How¬ 
ever, the photon spectrum measured by the WA98 Col¬ 
laboration showed a significant yield of direct photons at 
low transverse momentum. Thus the scenario 2 can be 
ruled out. Consequently, one has to assume that the ob¬ 
served direct photons of low pp had a significant elliptic 
anisotropy V 2 - of the same order of magnitude as the 
hadronic flow. The interpretation [4l|, of the low-pp 
direct photon yield measured by WA98 as dominantly 
produced by the bremsstrahlung process in the mesonic 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 


p.^ [GeV/c] 


FIG. 19: (Color on-line) Elliptic flow V 2 versus transverse mo¬ 
mentum Pt for the inclusive photons produced in 0-40% cen¬ 
tral Pb-|-Pb collisions at y^sivjv = 2.76 TeV from the PHSD 
(solid red line) in comparison to the data from the ALICE Col¬ 
laboration [331 ]: the blue error band reflects the finite statistics 
and the uncertainty in the modeling of the cross sections for 
the individual channels. 


collisions 7 r-|- 7 r—i-Tr-l-Tr-l-Tisin accord with the WA98 
data on the inclusive photon ^ 2 ( 7 *"'^^). 

Let us note that the same conclusions apply also to the 
most recent studies of the photon elliptic flow at RHIC 
and LHC. The PHENIX and ALICE Collaborations have 
measured the inclusive photon V 2 and found that at low 
transverse momenta it is comparable to the V 2 {pt) of de¬ 
cay photons as calculated in cocktail simulations based 
on the known mesonic V 2 {pt)- Therefore (a) either the 
yield of the direct photons to the inclusive ones is not 
statistically significant in comparison to the decay pho¬ 
tons or (b) the elliptic flow of the direct photons must be 
as large as ^ 2 ( 7 '^““^) and U 2 ( 7 ™'^ 0 - 

In Refs. 0 we have shown the elliptic flow of the 
inclusive and direct photons produced in minimum bias 
Au-f Au collisions at ^/snn = 200 GeV from the PHSD 
in comparison to the data of the PHENIX Collabora¬ 
tion. We found that the data on the inclusive photon 
V 2 could approximately be described. Furthermore, the 
pion decay photons dominate the inclusive photon spec¬ 
trum. Since the elliptic flow of pions is under control in 
PHSD in comparison to the data from th e PHENIX and 
STAR Collaborations (cf. Refs. [H. 1^. li^ . llOSL IIOGI ]). the 
spectrum of decay photons is also predicted reliably by 
the model. However, the good agreement with the in¬ 
clusive photon spectrum is especially meaningful due to 
the good description of the direct photon spectrum which 
was presented above in section [TV] 

After subtracting the contribution of the decay pho¬ 
tons from the inclusive photons, the direct photon V 2 
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FIG. 20: (Color on-line) Elliptic flow V 2 versus transverse mo¬ 
mentum pt for the direct photons produced in 0-40% central 
Pb-|-Pb collisions at ^Snn = 2.76 TeV from the PHSD (solid 
red line) in comparison to the preliminary data from the AL¬ 
ICE Collaboration [^; the blue error band is dominated by 
the uncertainty in the modeling of the cross sections for the 
individual channels. 


is accessed experimentally. In the PHSD, we calculate 
the direct photon ^ 2 ( 7 '^*’’) by building the weighted sum 
of the channels, which are not subtracted by the data- 
driven methods, as follows: the photons from the quark- 
gluon plasma, from the initial hard parton collisions 
(pQCD photons), from the decays of short-living reso¬ 
nances (oi-meson, (/)-meson, A-baryon), from the binary 
meson-l-meson and meson-l-baryon channels (tt -|- p —>■ 
TT-fy, 7r-|-7r — >■ p-l-7, V -\- p/n —>■ n/p -\- 7), and 
from the bremsstrahlung in the elastic meson-fmeson 
and meson-bbaryon collisions (to -|-to —> TO-I-TO-I-7, 
to-I-H—7 >to-|-H-|-7 ). We calculate the direct photon 
by summing up the elliptic flow of the individual chan¬ 
nels contributing to the direct photons, using their con¬ 
tributions to the spectrum as the relative pr-dependent 
weights, Wi{pt), i.e. 


'^ 2 ( 7 *’') = ^V 2 {f)w^{pT) 

i 


T,iM'l’‘)NiiPT) 


The results for the elliptic flow V 2 {pt) of direct pho¬ 
tons produced in Au + Au collisions at the top RHIC 
energy are shown in Fig. 1161 In comparison to the previ¬ 
ous results within the PHSD approach ia , the elliptic 
flow in the intermediate region of the transverse momenta 
1.0 < Pt < 2.0 GeV is reduced by about 50% due to the 
modifications in our computation of the bremsstrahlung 
channels beyond the soft-photon approximation in the 
present work. According to our calculations of the direct 
photon spectra (as presented above), almost half of the 


direct photons measured by PHENIX (in central colli¬ 
sions) stems from the collisions of quarks and gluons in 
the deconfined medium created in the initial phase of the 
collision. The photons produced in the QGP carry a very 
small V 2 and lead to an overall direct photon V 2 about a 
factor of 2 below the pion 'i;2(7r) even though the other 
channels in the sum (H5)l have large elliptic flow coeffi¬ 
cients V 2 of the order of ^ 2 ( 71 ) (cf. Fig. 7 of Ref. 0 )- 
Indeed, the parton collisions - producing photons in 
the QGP - take place throughout the evolution of the 
collision but the collision rate falls rapidly with time and 
thus the production of photons from the QGP is dom¬ 
inated by the early times. We illustrate this in FielTTl 
where it is seen that the rate of photon emission dN^/dt 
from QGP sources drops by orders of magnitude after 
the first 10 fm/c. As a consequence, the elliptic flow 
‘picked up’ by the photons from the parent parton col¬ 
lisions saturates after about a few fm/c and reaches a 
relatively low value of about 0.02, only. We note that 
a delayed prod uction of charges from the strong gluon 
fields (‘glasma’ 108l - llll |) might shift the QGP photon 
production to somewhat later times when the elliptic flow 
is built up more. However, we cannot quantitatively ans¬ 
wer whether the additional evolution in the pre-plasma 
state could generate considerable additional direct pho¬ 
ton V2- 


In the following we present calculations/predictions for 
the elliptic flow of inclusive and direct photons produced 
in Pb + Pb collisions at the energy of ^/snn = 2.76 TeV. 
First, we show in Fig. [18] the elliptic flow V 2 versus trans¬ 
verse momentum px for the charged hadrons produced in 
20-30% central Pb-|-Pb collisions at ^/snn = 2.76 TeV 
from the PHSD (solid red line); the blue band reflects 
the statistical uncertainty of the calculations. The PHSD 
results are c ompared to the data of the ALIGE Gollab- 
oration [Iq3 which suggests that the bulk dynamics is 
reasonably under control in the PHSD at these energies 
and the elliptic flow of the final charged particles is de¬ 
scribed up to Pt ~ 2 GeV/c in this centrality range. 

We, furthermore, present our calculations/predictions 
for the elliptic flow of inclusive and direct photons pro¬ 
duced in Pb + Pb collisions at the energy of ^/snn = 
2.76 TeV at the LHC within the acceptance of the AL¬ 
ICE detector. Fig. |T9| shows PHSD calculations for the 
elliptic flow V 2 versus transverse momentum px for the 
inclusive photons produced in 0-40% central Pb+Pb col¬ 
lisions at y/SNN = 2.76 TeV (solid red line) with the blue 
error band reflecting the finite statistics and the theoret¬ 
ical uncertainty in the modeling of the cross sections. A 
comparison to the respective data from the ALICE Col¬ 
laboration shows a comparable agreement as in case 
of the V 2 for the charged hadrons. 

The elliptic flow V 2 (px) of direct photons produced in 
0-40% central Pb-|-Pb collisions at ^/snn = 2.76 TeV 
from the PHSD (solid red line) is shown in Fig. |20| in 
comparison to the preliminary data from the ALICE Col¬ 
laboration [ 3 ^ : the blue error band is again dominated 
by the uncertainty in the modeling of the cross sections 
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FIG. 21: Centrality dependence of the direct photon V 2 f or Au +Au collisions at ^/sNN = 200 GeV for different centralities 
(see legend); the data from the PHENIX Gollaboration f^.[i07| are compared to the earlier PHSD predictions from Ref. 0]. 



Pt |GeV/c] 


FIG. 22: (Color on-line) Elliptic flow V 2 of direct photons from 
PHSD versus transverse momentum px produced in Pb-|-Pb 
collisions at -^snn = 2.76 TeV for different centrality classes 
(see legend). 


for the individual channels. As in case of Au-|-Au colli¬ 
sions at the top RHIC energy we slightly underestimate 
the direct photon V 2 , however, new data with higher ac¬ 
curacy will be needed to shed further light on the direct 
photon V 2 puzzle. 

The centrality dependence of the direct photon ellip¬ 
tic flow has been calculated within the PHSD in Ref. [3] 


for Au-|-Au collisions at the top RHIC energy. Fig. 
presents a direct comparison of the PHSD predictions for 
V 2 {pt) from Ref. 0] for Au-|-Au collisions at y/sNN=200 
GeV in the centrality classes 0-20% (a)^0-40% (b) and 
40-60% (c) with the data from Refs. Whereas 

the elliptic flow is roughly described in the most central 
class there is an increasing tendency to underestimate in 
the PHSD the strong elliptic flow especially for peripheral 
collisions where some additional source might be present. 
Thus the observed centrality dependence of the elliptic 
flow is roughly in agreement with the interpretation that 
a large fraction of the direct photons is of hadronic origin 
(in particular from the bremsstrahlung in meson-l-meson 
and meson-bbaryon collisions); the latter contribution be¬ 
comes stronger in more peripheral collisions. But more 
precise data will be mandatory for a robust conclusion. 

We close this Subsection in providing predictions for 
the centrality dependence of the direct photon V 2 {pt) in 
Pb-I-Pb collisions at ^/snn = 2.76 TeV in the centrality 
classes 0-20%, 20-40% and 40-80% which are of relevance 
for the upcoming measurements by the ALICE Collab¬ 
oration at the LHC. The actual results from PHSD are 
displayed in Fig. [52] and show a very similar centrality 
dependence as in case of Au-|-Au collisions at the top 
RHIC energy. 


C. Triangular flow of direct photons 

We have seen in the previous sections that the mea¬ 
sured spectra of direct photons could be reproduced by 
the PHSD calculations within a factor of 2 (which is com¬ 
parable with the current accuracy of the measurements). 
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FIG. 23: (Color on-line) Triangular flow V 3 versus transverse momentum pT for the direct photons produced in Au+Au collisions 
at ^^/snn = 200 GeV in three centrality classes (see legends). The PHSD res ults are s hown by the solid red lines in comparison 
to the data of the PHENIX Collaboration (black symbols) taken from Ref. [ioa,[n3. 


Also, the inclusive photon V 2 was well described and the 
elliptic flow of direct photons was qualitatively in line 
with the data (within a factor of 2 ) and attributed essen¬ 
tially to hadronic sources. 

On the other hand, there exists an alternative inter¬ 
pretation of the strong elliptic flow of direct photons, in 
which the azimuthal asymmetry of the photons is due 
to the initial strong magnetic field essentially produced 
by spectator charges (protons). Indeed, the magnetic 
field strength in the very early reaction stage reaches up 
to eBy K. hrn^ in semi-peripheral Au -\- Au collisions at 
^snn = 200 GeV (see the calculations within the PHSD 
in Ref. [H^; comparab le estimates have been obtained 
also in Refs. [mini)- These strong magnetic fields 
might influence the photon production via the polariza¬ 
tion of the medium, e.g. by influencing the motion of 
charged quarks in the QGP, or by directly inducing a real 
photon radiation via the virtual photon (R-field) coup¬ 
ling to a quark loop and (multiple) gluons; the photons 
are then produced azimuthal asymmetrically. 

The photon production under the influence of st rong 
mag netic fields has been calculated in Refs. [13, Ill4 
Ill6j| . The observed spectra and elliptic flow of direct 
photons could be explained using suitable assumptions 
on the conductivity, bulk viscosity or degree of chemical 
equilibration in the early produced matter. The common 
feature of these calculations was that the triangular flow 
coefficient W 3 of the direct photons was expected to be 
very small. Indeed, the magnetic field may lead to an 
azimuthal asymmetry V 2 but not to a triangular mode. 

Consequently, it is of interest to measure experimen¬ 
tally the third flow coefficient V 3 {pt) and to compare it 
to the calculations in the different classes of models: (a) 
those attributing the large elliptic flow and strong yield 



FIG. 24: (Color on-line) Triangular flow V 3 versus transverse 
momentum pr for the direct photons produced in different 
centrality classes for Pb-|-Pb collisions at ^ysNN = 2.76 TeV 
from the PHSD (see legend); the blue band reflects the uncer¬ 
tainty in the modeling of the cross sections for the individual 
channels and give a measure of the present level of accuracy. 


of direct photons to hadronic sources, e.g. the PHSD 
transport approach; (b) the models suggesting the large 
azimuthal asymmetry and additional yield of direct pho¬ 
tons to be caused by the early magnetic fields; (c) the 
models assuming that the yield of direct photons at low 
Pt is dominated by partonic channels. 

In Fig. [231 we present our results for the triangular flow 
U 3 versus transverse momentum pp for the direct pho¬ 
tons produced in Au-|-Au collisions at ^ysNN = 200 GeV 
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from the PHSD (solid red lines) for 0-20% (a), 20-40% 
(b) and 40-60% (c) centrality. The PHSD gives a posi¬ 
tive non-zero triangular flow of direct photons up to 6 % 
with very little centrality dependence on the level of the 
present accuracy 25%). The PHSD results are in 
agreement with the data of the PHENIX Collaboration 
from Refs. |ll2l . IllTl IllSjl which suggests that the sce¬ 
nario (a) is compatible with the measurements. 

Accordingly, the data from PHENIX in Eig. [53] do 
not point towards an interpretation of the direct photons 
being dominantly produced in the early stage under the 
influence of the magnetic field (b), because the of these 
photons is expected to be close to zero. Of course, the 
photon production in the magnetic fields occurs on top 
of other channels, which may carry finite V 3 . But the 
weighted sum of all the channels including the magnetic- 
field-induced photons will give a smaller U3 ^ 0 than 
the sum without this channel. The scenario (c) has been 
studied by other groups withi n a h ydrodynamic modeling 
of the collision in Refs. [n [m. The t riang ular flow 
V 3 {pt) of direct photons from Refs, [l^ iliof is about 
a factor of 2 smaller than that obtained in the present 
work from the PHSD. It will be possible to differentiate 
between the scenarios in the future, when data of higher 
accuracy will become available. 

We finally show in Pig. |53|our predictions for the trian¬ 
gular flow V 3 {pt) of direct photons produced in Pb + Pb 
collisions at ^/snn = 2.76 TeV for different centrali¬ 
ties as will be measured at the LHC. We expect the 
triangular flow of direct photons to reach at maximum 
ymax _ Q Q being comparable to the 

charged hadron ^3 [^. The centrality dependence of 
v^^pt) turns out to be low and is practically constant 
within the accuracy of the PHSD calculations. An ex¬ 
perimental confirmation of this expectation could further 
affirm the notion of large hadronic contributions to the 
direct photons and in particular the photon production 
via the bremsstrahlung in meson and baryon collisions. 


V. SUMMARY 

In the present work we have calculated the trans¬ 
verse momentum spectra, the elliptic flow vi and tri¬ 
angular flow V 3 of direct photons produced in Au-|-Au 
collisions at y/SNN = 200 GeV and in Pb + Pb colli¬ 
sions at ^/sNN = 2.76 TeV using the microscopic PHSD 
transport approach. Eor thermal photon production we 
have considered the interactions of quarks and gluons 
in the strongly interacting quark-gluon plasma (sQGP) 
(g + 9 -t 5 + 7 and q{q) + g ^ q{q) + 7 ), the pho¬ 
ton production in the hadronic decays (tt —>■ 7 -|- 7, 

5 —5>7-l-7, w—>'7r-|-7, 77 '—)>p-|-7, (/)—)>?7-|-7, 

ai —?> TT -I- 7, A —>■ V -f 7) as well as the interactions 
(TT-t-TT—>• p- 1 - 7 , V + TT ^ n + ^, the 2 — >• 2 interaction 
of mesons and baryons p + n/p ^ j + n/p, and last but 
not least the bremsstrahlung radiation in meson-f meson 
and meson-f baryon scattering m + m/ B —>■ m-\-m/B + ^ 


throughout the evolution of the collision. The pQGD 
photons produced in the initial hard binary scatterings 
are added to the thermal photons in order to obtain the 
final direct photon spectrum. 

In extension of our previous works on photon produc¬ 
tion in heavy-ion collisions at the top RHIC energy we 

(i) went beyond the soft-photon approximation (SPA) in 
the calculation of the bremsstrahlung processes meson + 
meson —>■ meson + meson + 7 , meson + baryon —>■ 
meson + baryon + 7 , 

(ii) quantified the suppression at low pr due the Landau- 
Pomeranchuk-Migdal (LPM) effect in connection with 
the electric conductivity, and 

(iii) incorporated the V -f A —>■ 7 -f A and A —> A -f 7 
channels into the PHSD. Furthermore, 

(iv) have presented calculations for Pb-fPb collisions at 
SPS and LHC energies. 

The result from Ref. [1^ for the reaction tt + tt ^ 
TT-f TT-f 7 , obtained within the one-boson exchange model 
beyond the soft-photon approximation - available up to 
the photon energy of 0.4 GeV - is confirmed by our 
present calculations that extend to 2 GeV. Furthermore, 
the improved SPA (1551) gives a very good approxima¬ 
tion to the exact result even at high y/s as long as the 
model for the elastic cross section is sufficiently realistic. 
In comparison, the constant-cross-section approximation 
based on the formula ((501) (used before) overestimates 
the exact rates for go > 1 GeV and underestimates for 
go < 0.4 GeV. The OBE differential cross sections have 
been implemented in the PHSD transport approach for 
photon production studies in heavy-ion reactions from 
SPS to LHC energies. 

Furthermore, we have given an estimate for the pho¬ 
ton suppression due to the Landau-Pomeranchuk-Migdal 
(LPM) effect employing the approximation (l 36 ll that cap¬ 
tures the relevant physics. When inserting the interac¬ 
tion time of the partons from PHSD at finite temperature 
we obtain a photon rate that is well in line with Eq. dSIl) 
in the limit of vanishing photon energy. We recall that 
the electric conductivity uq regulates the divergence of 
the photon rate and has been evaluated in the PHSD in 
Ref. ffT\ in rough agreement with related results from 
lattice QCD. The actual results in Fig. 9 demonstrate 
that the LPM effect becomes important for photon ener¬ 
gies below 0.4 GeV in case of dense and strongly inter¬ 
acting partons. Since the experimental photon spectra 
measured so far at RHIC and LHC energies start from 
higher phton energies the LPM effect might be safely dis¬ 
carded. 

In case of relativistic heavy-ion collisions we have found 
that the PHSD calculations reproduce the transverse mo¬ 
mentum spectra of direct photons as measured by the 
PHENIX Collaboration in Refs. [ 2 ^, [s^ for Au-|-Au col¬ 
lisions at the top RHIC energy. The calculations reveal 
the channel decomposition of the observed direct photon 
spectrum and show that the photons produced in the 
QGP constitute at most about 50% of the direct photons 
in central collisions with the rest being distributed among 
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the other channels: mesonic interactions, decays of mas¬ 
sive hadronic resonances and the initial hard scatterings. 
Our calculations demonstrate that the photon produc¬ 
tion in the QGP is dominated by the early phase (similar 
to hydrodynamic models) and is localized in the center 
of the fireball, where the collective flow is still rather low, 
i.e. on the 2-3 % level, only. Thus, the strong V 2 of di¬ 
rect photons - which is comparable to the hadronic V 2 - in 
PHSD is attributed to hadronic channels, i.e. to meson 
and baryon induced reactions. On the other hand, the 
strong V 2 of the ’parent’ hadrons, in turn, stems from the 
interactions in the QGP via collisions and the partonic 
mean-field potentials. Accordingly, the presence of the 
QGP shows up ’indirectly’ in the direct photon elliptic 
flow. 

We, furthermore, have demonstrated that the ellip¬ 
tic flow of charged hadrons from PHSD is in a rea¬ 
sonable agreement with the ALICE data. Predic¬ 
tions/calculations for the inclusive photon V 2 and direct 
photon V 2 have been provided as well as associated re¬ 
sults for the triangular direct photon flow v^{pt) which 
is as large as the charged hadron at = 200 GeV 

and y/SNN = 2.76 TeV. The large triangular flow of di¬ 
rect photons seen by PHENIX at the top RHIC energy is 
in line with the PHSD calculations and in conflict with 
a scenario that attributes the large direct photon V 2 to 
the influence of the early strong magnetic field from the 
spectator charges. 

The centrality dependence of the direct photon yield 
and flow has the potential to further clarify the direct 
photon production mechanisms. We find a good agree¬ 
ment between the PHENIX measurements and the PHSD 
calculations at the top RHIC energy. In particular, the 
integrated thermal photon yield in PHSD was predicted 
to scale as with the exponent a ~ 1.5, which is in 

a good agreement with the most recent measurement of 


a = 1.48±0.08±0.04 by the PHENIX Collaboration [HI. 
This observation supports the conclusion that the low 
transverse momentum direct photons have a strong con¬ 
tribution from the binary hadronic photon production 
sources, such as the meson + meson and meson + baryon 
bremsstrahlung. It will be important to investigate ex¬ 
perimentally the scaling of the direct photon yield and 
flows V 2 and V 3 with centrality also in Pb-|-Pb collisions 
at the LHC energies for which we have provided explicit 
predictions from the PHSD. 
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